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Abstract 


These lecture notes provide a self-contained introduction to the mathematical methods required in a Bachelor 
degree programme in Business, Economics, or Management. In particular, the topics covered comprise real¬ 
valued vector and matrix algebra, systems of linear algebraic equations, Leontief’s stationary input-output 
matrix model, linear programming, elementary financial mathematics, as well as differential and integral 
calculus of real-valued functions of one real variable. A special focus is set on applications in quantitative 
economical modelling. 
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Qualification objectives of the module 
(excerpt) 


The qualification objectives shall be reached by an integrative approach. 

A broad instructive range is aspired. The students shall acquire a 360 degree orientation concerning 
the task- and personnel-related tasks and roles of a manager and supporting tools and methods and 
be able to describe the coherence in an integrative way. The knowledge concerning the tasks and 
the understanding of methods and tools shall be strengthened by a constructivist approach based 
on case studies and exercises. 

Students who have successfully participated in this module will be able to 


• solve problems in Linear Algebra and Analysis and apply such mathematical methods to 
quantitative problems in management. 

• apply and challenge the knowledge critically on current issues and selected case studies. 
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Introduction 


These lecture notes contain the entire material of the quantitative methods part of the first semester 
module 0.1.1 IMQM: Introduction to Management and its Quantitative Methods at Karl- 
shochschule International University. The aim is to provide a selection of tried-and-tested math¬ 
ematical tools that proved efficient in actual practical problems of Economics and Management. 
These tools constitute the foundation for a systematic treatment of the typical kinds of quantita¬ 
tive problems one is confronted with in a Bachelor degree programme. Nevertheless, they provide 
a sufficient amount of points of contact with a quantitatively oriented subsequent Master degree 
programme in Economics, Management, or the Social Sciences. 

The prerequisites for a proper understanding of these lecture notes are modest, as they do not go 
much beyond the basic A-levels standards in Mathematics. Besides the four fundamental arith¬ 
metical operations of addition, subtraction, multiplication and division of real numbers, you should 
be familiar, e.g., with manipulating fractions, dealing with powers of real numbers, the binomial 
formulae, determining the point of intersection for two straight lines in the Euclidian plane, solv¬ 
ing a quadratic algebraic equation, and the rules of differentiation of real-valued functions of one 
variable. 

It might be useful for the reader to have available a modem graphic display calculator (GDC) for 
dealing with some of the calculations that necessarily arise along the way, when confronted with 
specific quantitative problems. Some current models used in public schools and in undergraduate 
studies are, amongst others, 

• Texas Instruments TI-84 plus, 

• Casio CFX-9850GB PLUS. 

However, the reader is strongly encouraged to think about resorting, as an alternative, to a spread¬ 
sheet programme such as EXCEL or OpenOffice to handle the calculations one encounters in 
one’s quantitative work. 

The central theme of these lecture notes is the acquisition and application of a number of effec¬ 
tive mathematical methods in a business oriented environment. In particular, we hereby focus on 

quantitative processes of the sort 


INPUT ^ OUTPUT , 

for which different kinds of functional relationships between some numerical INPUT quantities 
and some numerical OUTPUT quantities are being considered. Of special interest in this context 
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will be ratios of the structure 

OUTPUT 
INPUT ■ 

In this respect, it is a general objective in Economics to look for ways to optimise the value of 
such ratios (in favour of some economic agent), either by seeking to increase the OUTPUT when 
the INPUT is confined to be fixed, or by seeking to decrease the INPUT when the OUTPUT is 
confined to be fixed. Consequently, most of the subsequent considerations in these lecture notes 
will therefore deal with issues of optimisation of given functional relationships between some 
variables, which manifest themselves either in minimisation or in maximisation procedures. 

The structure of these lecture notes is the following. Part I presents selected mathematical methods 
from Linear Algebra, which are discussed in Chs.[I]to[5l Applications of these methods focus on 
the quantitative aspects of flows of goods in simple economic models, as well as on problems in 
linear programming. In Part II, which is limited to Ch. we turn to discuss elementary aspects 
of Financial Mathematics. Fundamental principles of Analysis, comprising differential and in¬ 
tegral calculus for real-valued functions of one real variable, and their application to quantitative 
economic problems, are reviewed in Part III; this extends across Chs.|7]and[8l 

We emphasise the fact that there are no explicit examples nor exercises included in these lecture 
notes. These are reserved exclusively for the lectures given throughout term time. 

Recommended textbooks accompanying the lectures are the works by Asano (2013) [|2l, Dowl¬ 
ing (2009) [fTll . Dowling (1990) ifT^ . Bauer et al (2008) S, Bosch (2003) |l6l, and Hiilsmann 
et al (2005) ifT^ . Some standard references of Applied Mathematics are, e.g., Bronstein et al 
(2005) ATI and Arens et al (2008) [HI. Should the reader feel inspired by the aesthetics, beauty, 
ellegance and efficiency of the mathematical methods presented, and, hence, would like to know 
more about their background and relevance, as well as being introduced to further mathematical 
techniques of interest, she/he is recommended to take a look at the brilliant books by Penrose 
(2004) fflU . Singh (1997) [l23l . Gleick(1987) [IT^ and Smith (2007) [l24l . Note that most of the 
textbooks and monographs mentioned in this Introduction are available from the library at Karl- 
shochschule International University. 

Finally, we draw the reader’s attention to the fact that the *.pdf version of these lecture notes con¬ 
tains interactive features such as fully hyperlinked references to original publications at the web¬ 
sites dx . doi . org and j stor . org, as well as active links to biographical information on sci¬ 
entists that have been influential in the historical development of Mathematics, hosted by the web¬ 
sites The MacTutor History of Mathematics archive (www-hi story . mcs . st-and. ac . uk) 
and en . wikipedia . org, 




Chapter 1 

Vector algebra in Euclidian space 


Let us begin our elementary considerations of vector algebra with the introduction of a special 
class of mathematical objects. These will be useful at a later stage, when we turn to formulate 
certain problems of a quantitative nature in a compact and elegant way. Besides introducing these 
mathematical objects, we also need to define which kinds of mathematical operations they can be 
subjected to, and what computational rules we have to take care of. 


1.1 Basic concepts 


Given be a set V of mathematical objects a which, for now, we want to consider merely as a 
collection of n arbitrary real numbers oi,..., a*,..., a^. In explicit terms. 



< 

/ ai \ 


> 

V = < 

a = 

(Xi 

Oj G M, f = 1,.. 

., n > 


< 

V j 


> 


Formally the n real numbers considered can either be assembled in an ordered pattern as a column 
or a row. We define 

Def.r Real-valued column vector with n components 



( 1 . 2 ) 


nxl 


Notation: a G M’ 
and 
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CHAPTER 1. VECTOR AEGEBRA IN EUCEIDIAN SPACE 
Def.: Real-valued row vector with n components 


a 


(oi, 


cti G 


i = 1, 


I n 


(1.3) 


Notation: G 

Correspondingly, we define the n-component objects 


/o\ 


0 ;= 


and 0^ := (0,..., 0,..., 0) 


Voy 


(1.4) 


to constitute related zero vectors. 

Next we define for like objects in the set V, i.e., either for n-component column vectors or for 
n-component row vectors, two simple computational operations. These are 

Def.: Addition of vectors 



(1.5) 


and 

Def.: Rescaling of vectors 



( 1 . 6 ) 


The rescaling of a vector a with an arbitrary non-zero real number A has the following effects: 

• IAI > 1 — stretching of the length of a 

• 0 < IAI < 1 — shrinking of the length of a 

• A < 0 — directional reversal of a. 
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1.2. DIMENSION AND BASIS OF 

The notion of the length of a veetor a will be made preeise shortly. 

The addition and the rescaling of n-component vectors satisfy the following addition and multipli¬ 
cation laws: 

Computational rules for addition and rescaling of vectors 

For vectors a, b,c G M”: 

1. a + b = b + a (commutative addition) 

2. a + {b + c) = {a + b) + c (associative addition) 

3. a + 0 = a (addition identity element) 

4. For every a,b G M”, there exists exactly one a; G such that a + x = b 

(invertibility of addition) 

5. (A/i)a = X{fia) with A G M (associative rescaling) 

6. la = a (rescaling identity element) 

7. A(ci T b) = Xa T Abj 

(A + fi)a = Xa + fia with A, /r G M (distributive rescaling). 

In conclusion of this section, we remark that every set of mathematical objects V constructed in 
line with Eq. (11.11) . with an addition and a rescaling defined according to Eqs. (11.51) and (11.61). and 
satisfying the laws stated above, constitutes a linear vector space over Euclidian space 

1.2 Dimension and basis of 

Let there be given m n-component vector^ ai,..., a,,..., am G M”, as well as m real numbers 
Ai,..., Aj,..., Am G M. The new n-component vector b resulting from the addition of arbitrarily 
rescaled versions of these m vectors according to 


m 

b = Aid-i -f ... -f Xj^a^ -f ... -f Xmf^m ^ ^ AjUj G 

i=l 


(1.7) 


is referred to as a linear combination of the m vectors ai, i = 1,. ■ ■ ,m. 

Def.: A set of m vectors ai,..., a,,..., am G M" is called linearly independent when the condi¬ 
tion 

m 

0 = Aitti -f . . . -f AiUj -f . . . -f Xmdm = ^ ^ AiUj , (1-8) 

i=l 

'This is named after the ancient greek mathematician Euclid of Alexandria (about 325 BC-265 BC) 

slightly s horter notation for n-component column vectors a G is given by a G R"; likewise G R" 

for n-component row vectors G 
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i.e., the problem of forming the zero vector 0 G M” from a linear eombination of the m veetors 
ai,..., Oj,..., £ IK”, can only be solved trivially, namely by 0 = Ai = ... = Aj = ... = A^. 

When, however, this eondition ean be solved non-trivially, with some A* ^ 0, then the set of m 
veetors ai,..., a*,..., G is ealled linearly dependent. 

In Euelidian spaee M”, there is a maximum number n (!) of veetors whieh ean be linearly indepen¬ 
dent. This maximum number is referred to as the dimension of Euclidian space M”. Every set 
of n linearly independent veetors in Euelidian spaee M” eonstitutes a possible basis of Euclidian 
space M”. If the set {oi, a„} eonstitutes a basis of M”, then every other veetor b G M” 

ean be expressed in terms of these basis veetors by 


b — + . . . -f /Sitti + . . . -f /SnCLn — • (1-9) 

i=l 

The resealing faetors /S* G M of the a, G M” are ealled the components of vector b with respect 
to the basis {ai, ..., a*, ..., a„}. 


Remark: The n unit vectors 


(^ \ 


/o\ 


/o\ 

0 

, 62 := 

1 

, . . . , 6n . 

0 

■ 0 


■ 0 






( 1 . 10 ) 


constitute the so-called canonical basis of Euclidian space M”. With respect to this basis, all 
vectors 6 G can be represented as a linear combinationen 


b = 




= biEi + 6262 H-h 


\bn J 


n 

i=\ 


( 1 . 11 ) 


1.3 Euclidian scalar product 

Einally, to conclude this section, we introduce a third mathematical operation defined for vectors 
on R”. 

Def.: Eor an n-component row vector G R^^” and an n-component column vector b G 

the Euclidian scalar product 



( 1 . 12 ) 




1.3. EUCLIDIAN SCALAR PRODUCT 
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defines a mapping )-R from the produet set of n-eomponent row and eolumn 

veetors to the set of real numbers. Note that, in eontrast to the addition and the resealing of n- 
eomponent veetors, the outeome of forming a Euelidian sealar produet between two n-eomponent 
veetors is a single real number. 

In the eontext of the Euelidian sealar produet, two non-zero veetors a, 6 G M” (wit a ^ 0 ^ b) 
are referred to as mutually orthogonal when they exhibit the property that 0 = ■ b = ■ a. 

Computational rules for Euclidian scalar product of vectors 

Eor vectors a,b,c G M”: 

1. {a + 6)^ ■ c = ■ c + b^ ■ c (distributive scalar product) 

2. ■ b = b^ ■ a (commutative scalar product) 

3. (Aa^) ■ b = X{a^ ■ b) with A G M (homogeneous scalar product) 

4. a > 0 for all a 0 (positive definite scalar product). 

Now we turn to introduce the notion of the length of an n-component vector. 

Def.: The length of a vector a G M” is defined via the Euclidian scalar product as 



n 

\a\ := ■ a = Jal + ... + af + ... + al =■ 


\ 

i=l 


(1.13) 


Technically one refers to the non-negative real number |a| as the absolute value or the Euclidian 
norm of the vector a G M”. The length of a G M” has the following properties: 

• |a| > 0, and |a| = 0 -v^ a = 0; 

• |Aa| = |A||a| for A G M; 

• |a-l-6|<|a|-|-|b| (triangle inequality). 

Every non-zero vector a G M”, i.e., |a| > 0, can be rescaled by the reciprocal of its length. This 
procedure defines the 

Def.: Normalisation of a vector a G M”; 

a := ^ |d,| = 1 . (1-14) 

a 


By this method one generates a vector of length 1, i.e., a unit vector a. To denote unit vectors we 
will employ the “hat” symbol. 

Easily, also by means of the Euclidian scalar product, we introduce the angle enclosed between 
two non-zero vectors. 
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Def.: Angle enclosed between a,b ^ 0 & R^ 


cos[(^(a, b)] 



(^(a, b) = cos ■ b) . 


(1.15) 


Remark: The inverse cosine functiorEI cos ^(...) is available on every standard GDC or spread¬ 
sheet. 


^The notion of on inverse function will be discussed later in Ch.|7] 






Chapter 2 
Matrices 


In this chapter, we introduce a second class of mathematical objects that are more general than 
vectors. For these objects, we will also define certain mathematical operations, and a set of com¬ 
putational rules that apply in this context. 


2.1 Matrices as linear mappings 


Consider given a collection of m x n arbitrary real numbers an, au ■ ■ ■, a-ij, ■ ■ ■, a-mn, which we 
arrange systematically in a particular kind of array. 


Def.r A real-valued (m x n)-matrix is formally defined to constitute an array of real numbers 
according to 


an 

0-12 

aij 

^In 

0-21 

0-22 

... a2j ... 

Ct2n 

ttii 

Oi2 

. . . Ctij . . . 

^in 


Om2 

• • • ^mj • • • 

^mn 


where aij G M, f = 1,..., m; j = 1,..., n. 

Notation: A G 


( 2 . 1 ) 


Characteristic features of this array of real numbers are: 


• m denotes the number of rows of A, n the number of columns of A. 

• ttij represents the elements of A; aij is located at the point of intersection of the Ah row and 
the jth column of A. 

• elements of the Ah row constitute the row vector (an, ai 2 , • • •, aij, • • •, ciin), elements of the 
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jth column the column vector 


Formally column vectors need to be viewed as (n x 1) -matrices, row vectors as (1 x n) -matrices. 
An (m X n)-zero matrix, denoted by 0, has all its elements equal to zero, i.e., 

/ 0 0 ... 0 \ 

0 0 ... 0 

0 ;= ... . . ( 2 . 2 ) 
^0 0 ... 0 y 

Matrices which have an equal number of rows and columns, i.e. m = n, are referred to as 
quadratic matrices. In particular, the (n x n)-unit matrix (or identity matrix) 


/1 

0 . 

. 0 . 

• o\ 

• o 

1 . 

• o 

• o 

0 

0 . 

. 1 . 

. 0 

• • o 

0 . 

. 0 . 

. 1 / 


holds a special status in the family of (n x n) -matrices. 

Now we make explicit in what sense we will comprehend (m x n) -matrices as mathematical 
objects. 

Def.: A real-valued matrix A G defines by the computational operation 


Aa; : = 


ail ai2 
®21 ®22 

an ai2 


ami ®m2 • • • amj 

aiiXi -f ai2X2 + • • • 
a2lXi -\- 0222^2 + • • • 


ain \ / Xi 

a2n 1 I X2 


aijXj -\- 
a2jXj -f 


ainXn 

®2n2^n 


anXi -f ai2X2 -f ... -f aijXj -f ... -|- a^Xn 
amlO^l Um22^2 “f • • • “f amjXj -f . . . -f amnO^r 


= y (2.4) 
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a mapping A : —)■ i.e. a mapping from the set of real-valued n-eomponent eolumn 

veetors (here: x) to the set of real-valued m-eomponent eolumn veetors (here: y). 

In loose analogy to the photographic process, x can be viewed as representing an “object,” A a 
“camera,” and y the resultant “image.” 

Since for real-valued vectors Xi^X 2 G R"^^ and real numbers A G R, mappings defined by real¬ 
valued matrices A G R"^^" exhibit the two special properties 


A(a3i -f X 2 ) = (Acci) -f (Acca) 

A(Acci) = A(Aa;i) , 

they constitute linear mappingsQ 


(2.5) 


We now turn to discuss the most important mathematical operations defined for (m x n)-matrices, 
as well as the computational rules that obtain. 


2.2 Basic concepts 


Def.: Transpose of a matrix 

For A G R™^”, we define the process of transposing A by 


where i = 1,... ,m und j = 1,... ,n. Note that it holds that A^ G 


( 2 . 6 ) 


When transposing an (m x tt.) -matrix, one simply has to exchange the matrix’ rows with its columns 
(and vice versa): the elements of the first row become the elements of the first column, etc. It 
follows that, in particular, 

(A^)^ = A (2.7) 


applies. 


Two special cases may occur for quadratic matrices (where m = n): 


• When A^ = A, one refers to A as a symmetric matrix. 


• When A^ = —A, one refers to A as an antisymmetric matrix. 


Def.: Addition of matrices 

For A, B G R*"^”, the sum is given by 

A -f B =: C: ciij ~\~ h^j =: Cij , 


( 2 . 8 ) 


with z = 1,..., m and j = 1,..., n. 

Note that an addition can be performed meaningfully only for matrices of the same format. 

Tt is important to note at this point that many advanced mathematical models designed to describe quantitative 
aspects of some natural and economic phenomena do not satisfy the conditions (I2.51 l. as they employ non-linear 
mappings for this purpose. However, in such contexts, linear mappings often provide useful first approximations. 
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Def.: Rescaling of matrices 
For A G and A G M\{0}, let 


A A. —! C i \(lij —Cij , 


(2.9) 


where i = 1,...,m and j = 1,... ,n. 

When rescaling a matrix, all its elements simply have to be multiplied by the same non-zero real 
number A. 

Computational rules for addition and rescaling of matrices 

For matrices A, B, C G 

1. A + B = B-hA (commutative addition) 

2. A + (B + C) = (A B) + C (associative addition) 

3. A -h 0 = A (addition identity element) 

4. For every A and B, there exists exactly one Z such that A + Z = B. 


(invertibility of addition) 

5. (A/r)A = A(/rA) with A, /r G R\{0} (associative rescaling) 

6. 1A = A (rescaling identity element) 

7. A(A + B) = AA +AB; 

(A + f^)A = AA + nA with A, /r G R\{0} (distributive rescaling) 

8. (A B)^ = A^ + B^ (transposition rule 1) 

9. (AA)^ = AA^ with A G R\{0}. (transposition rule 2) 


Next we introduce a particularly useful mathematical operation for matrices. 


2.3 Matrix multiplication 

Def.: For a real-valued (m x n)-matrix A and a real-valued (n x r)-matrix B, a matrix multipli¬ 
cation is defined by 


AB =; C 

"F . . . “f (^ikbkj “F . . . “F (^in^nj . —i (^ik^kj . Qj ) 


( 2 . 10 ) 


with i = 1,... ,m and j = 1,... ,r, thus yielding as an outcome a real-valued (m x r)-matrix C. 
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The element of C at the interseetion of the ith row and the jth eolumn is determined by the 
eomputational rule 

Cij = Euelidian sealar produet of ith row veetor of A and jth eolumn veetor of B . (2.11) 

It is important to realise that the definition of a matrix multiplieation just provided depends in 
an essential way on the faet that matrix A on the left in the product needs to have as many (!) 
columns as matrix B on the right rows. Otherwise, a matrix multiplieation cannot be defined in a 
meaningful way. 

GDC: For matriees [A] and [B] edited beforehand, of matehing formats, their matrix multiplieation 
ean be evaluated in mode MATRIX —)■ NAMES by [A] * [B]. 

Computational rules for matrix multiplication 

For A, B, C real-valued matriees of eorrespondingly matehing formats we have: 


1. AB = 0 is possible with A 7 ^ 0, B 7 ^ 0. 

2. A(BC) = (AB)C 

3. A ^ ^ A = A 

gKnxn gjgmxm 

4. (A + B)C = AC + BC 
C(A + B) = CA + CB 

5. A(AB) = (AA)B = A(AB) with A G M 

6. (AB)^ = B^A'^ 


(zero divisor) 
(associative matrix multiplication) 
(multiplicative identity element) 

(distributive matrix multiplication) 
(homogeneous matrix multiplication) 
(transposition rule). 
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Chapter 3 

Systems of linear algebraic equations 


In this chapter, we turn to address a particular field of application of the notions of matrices and 
vectors, or of linear mappings in general. 


3.1 Basic concepts 

Let us begin with a system of m G N linear algebraic equations, wherein every single equation can 
be understood to constitute a constraint on the range of values of n G N variables xi,..., G M. 
The objective is to determine all possible values of Xi,..., G M which satisfy these constraints 
simultaneously. Problems of this kind, namely systems of linear algebraic equations, are often 
represented in the form 

• Representation 1: 


T • • • T T • . • T CLln^n — 

G>i\X\ T • • • “t“ ^ij^j T .. • T Oji'dX'fi hj^ (3.1) 

T . . • T O/rnjXj T • • • T * 

Depending on how the natural numbers m and n relate to one another, systems of linear algebraic 
equations can be classified as follows: 

• m < n: fewer equations than variables; the linear system is under-determined, 

• m = n: same number of equations as variables; the linear system is well-determined, 

• m > n: more equations than variables; the linear system is over-determined. 

A more compact representation of a linear system of format (m x n) is given by 
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• Representation 2: 



/ Oil ... ttij ... ain \ 


/ \ 


/ w \ 


Aa: = 

dii ... dij ... diYi 



= 

h 

= b. 


V • • • ^mj • • • ^mn ) 


\xn j 


\b^ ) 



The mathematical objects employed in this variant of a linear system are as follows: A takes the 
central role of the coefficient matrix of the linear system, of format {m x n), x is its variable 
vector, of format (n x 1), and, lastly, b is its image vector, of format (m x 1). 

When dealing with systems of linear algebraic equations in the form of Representation 2, i.e. 
Ax = b, the main question to be answered is: 

Question: For given coefficient matrix A and image vector b, can we find a variable vector x 
that A maps onto b? 

In a sense this describes the inversion of the photographic process we had previously referred to: 
we have given the camera and we already know the image, but we have yet to find a matching 
object. Remarkably, to address this issue, we can fall back on a simple algorithmic method due to 
the German mathematician and astronomer Carl Friedrich GauB (1777-1855). 


3.2 GauBian elimination 

The algorithmic solution technique developed by GauB is based on the insight that the solution set 
of a linear system of m algebraic equations for n real-valued variables, i.e. 



(3.3) 


remains unchanged under the following algebraic equivalence transformations of the linear sys¬ 
tem: 


1 . changing the order amongst the equations, 

2. multiplication of any equation by a non-zero real number c 7^ 0, 

3. addition of a multiple of one equation to another equation, 

4. changing the order amongst the equations. 


Specifically, this implies that we may manipulate a given linear system by means of these four 
different kinds of equivalence transformations without ever changing its identity. In concrete cases, 
however, one should not apply these equivalence transformations at random but rather follow a 
target oriented strategy. This is what GauBian elimination can provide. 
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Target: To cast the augmented coefficient matrix (A|6), i.e., the array 

Oil ... aij ... ain bi 

G>il • • • Giij . . . Ctin bi , (3.4-) 

Clrnl ■ ■ ■ Ojmj • • • ^mn 

when possible, into upper triangular form 


1 . 

. aij .. 

• • 5>in 

bi 


0 . 

Ctij . ■ 

• • ^in 

bi , 

(3.5) 

0 . 

. 0 

• • ^mn 

bm 



by means of the four kinds of equivalence transformations such that the resultant simpler final 
linear system may easily be solved using backward substitution. 

Three exclusive cases of possible solution content for a given system of linear algebraic equations 
do exist. The linear system may possess either 

1 . no solution at all, or 

2 . a unique solution, or 

3. multiple solutions. 

Remark; Linear systems that are under-determined, i.e., when m < n, can never be solved 
uniquely due to the fact that in such a case there not exist enough equations to constrain the values 
of all of the n variables. 

GDC: For a stored augmented coefficient matrix [A] of format (m x n +1), associated with a given 
(m X n) linear system, select mode MATRIX —)■ MATH and then call the function rref ([A]). It is 
possible that backward substitution needs to be employed to obtain the final solution. 

For completeness, we want to turn briefly to the issue of solvability of a system of linear algebraic 
equations. To this end, we need to introduce the notion of the rank of a matrix. 

3.3 Rank of a matrix 

Def.r A real-valued matrix A G possesses the rank 

rank(A) = r , r < min{m, n} (3.6) 

if and only if r is the maximum number of row resp. column vectors of A which are linearly 
independent. Clearly, r can only be as large as the smaller of the numbers m and n that determine 
the format of A. 
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For quadratic matrices A e there is available a more elegant measure to determine its 

rank. This (in the present ease real-valued) measure is referred to as the determinant of matrix A, 
det(A), and is defined as follows. 

Def.: 


(i) When A G its determinant is given by 


det(A) : = 


Oil ai2 
O 2 I ^22 


0-11022 ~ O12O21 , 


(3.7) 


i.e. the differenee between the produets of A’s on-diagonal elements and A’s off-diagonal 
elements. 


(ii) 


When A G the definition of A’s determinant is more eomplex. 

given by 



Oil 

fll2 

Ol3 

det(A) : = 

^21 

^22 

023 


^31 

032 

033 


In that ease it is 


011(022033 — O 32 O 23 ) + 021(032013 — O 12 O 33 ) -f 031(012023 — O 22 O 13 ) (.3.8) 


Observe, term by term, the eyclic permutation of the first index of the elements Ojj according 
to the rule 1 —)■ 2 —)■ 3 —)■ 1. 


(iii) Finally, for the (slightly involved) definition of the determinant of a higher-dimensional 
matrix A G please refer to the literature; e.g. Bronstein et al (2005) [[71 p 267]. 


To determine the rank of a given quadtratic matrix A G one now installs the following 

criteria: rank(A) = r = n, if det(A) 7 ^ 0, and rank(A) = r < n, if det(A) = 0. In the first case, 
A is referred to as regular, in the second as singular. For quadratic matrices A that are singular, 
rank( A) = r (with r < n) is given by the number r of rows (or columns) of the largest possible 
non-zero subdeterminant of A. 

GDC: For a stored quadratic matrix [A], select mode MATRIX —)■ MATH and obtain its determinant 
by calling the function det([A]). 


3.4 Criteria for solving systems of linear algebraic equations 

Making use of the concept of the rank of a real-valued matrix A G we can now summarise 

the solution content of a specific system of linear algebraic equations of format (m x n) in a table. 
For given linear system 

Ax = b , 

with coefficient matrix A G R'"^”, variable vector x G R"^^ and image vector b G R™^^, there 
exist(s) 
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by^O 

b = 0 

1. rank(A) ^ rank(A|b) 

no solution 


2. rank(A) = rank(A b) = r 



(a) r = n 

a unique 
solution 

X = 0 

(b) r < n 

multiple 
solutions: 
n — r free 
parameters 

multiple 
solutions: 
n — r free 
parameters 


(A|6) here denotes the augmented eoeffieient matrix. 

Next we diseuss a partieularly useful property of regular quadratie matriees. 


3.5 Inverse of a regular {n x n) -matrix 

Def.: Let a real-valued quadratie matrix A G be regular, i.e., det(A) G M\{0}. Then there 
exists an inverse matrix A“^ to A defined by the eharacterising properties 


A-^A = AA"^ = 1 . 


(3.9) 


Here 1 denotes the (n x n)-unit matrix [cf. Eq. (12.31) 1. 

When a eomputational deviee is not at hand, the inverse matrix A“^ of a regular quadratie matrix 
A ean be obtained by solving the matrix-valued linear system 

AX=1 (3.10) 

for the unknown matrix X by means of simultaneous GauBian elimination. 

GDC: For a stored quadratie matrix [A], its inverse matrix ean be simply obtained as [A]~\ where 
the x~^ funetion key needs to be used. 
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Computational rules for the inverse operation 

For A, B G with det(A) 7 ^ 0 7 ^ det(B), it holds that 

1. (A-i)-i = A 

2. (AB)-i = B-^A-i 

3. (A^)-i = (A-Y 

4. (AA)~i = 

A 

The special interest in applications in the concept of inverse matrices arises for the following 
reason. Consider given a well-determined linear system 

Acc = b , 

with regular quadratic coefficient matrix A G i.e., det(A) 7 ^ 0. Then, for A, there exists an 
inverse matrix A~^ Matrix-multiplying both sides of the equation above/rom the left (!) by the 
inverse A“\ results in 


A ^(Acc) = (A ^A)a; = la; = a:: = ■ (3-11) 

left-hand side right-hand side 

In this case, the unique solution (!) x = A~^b of the linear system arises simply from matrix 
multiplication of the image vector b by the inverse matrix of A. (Of course, it might actually 
require a bit of computational work to determine A~^.) 


3.6 Outlook 

There are a number of exciting advanced topics in Linear Algebra. Amongst them one finds 
the concept of the characteristic eigenvalues and associated eigenvectors of quadratic matrices, 
which has particularly high relevance in practical applications. The question to be answered here 
is the following: for given real-valued quadratic matrix A G R”^”, do there exist real numbers 
An G R and real-valued vectors Vn G R"'^^ which satisfy the condition 

AVn = XnVn^ (3-12) 

Put differently: for which vectors G R"^^^ does their mapping by a quadratic matrix A G R"^” 
amount to simple rescalings by real numbers An G R? 

By re-arranging, Eq. (13.121) can be recast into the form 


n 1 


0 = (A — Anl) V 


(3.13) 
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with 1 an (n X n)-unit matrix [cf. Eq. (12.31) 1 and 0 an n-component zero vector. This condition 
corresponds to a homogeneous system of linear algebraic equations of format {n x n). Non-trivial 
solutions 7 ^ 0 to this system exist provided that the so-called characteristic equation 

0 = det (A - A„l) , (3.14) 

a polynomial of degree n (cf. Sec. 17.1.11) . allows for real-valued roots A„ G M. Note that symmetric 
quadratic matrices (cf. Sec. 12.21) possess exclusively real-valued eigenvalues A„. When these 
eigenvalues turn out to be all different, then the associated eigenvectors Vn prove to be mutually 
orthogonal. 

Knowledge of the spectrum of eigenvalues An G M and associated eigenvectors Vn G 
of a real-valued matrix A G provides the basis of a transformation of A to its diagonal 
form Aa„, thus yielding a diagonal matrix which features the eigenvalues An as its on-diagonal 
elements; cf. Leon (2009) IfT^ . 

Amongst other examples, the concept of eigenvalues and eigenvectors of quadratic real-valued 
matrices plays a special role in Statistics, in the context of exploratory principal component 
analyses of multivariate data sets, where the objective is to identify dominant intrinsic structures; 
cf. Hair et al (2010) [|T4l Ch. 3] and Ref. IfT^ App. A]. 
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Chapter 4 

Leontief’s stationary input-output matrix 
model 


We now turn to discuss some specific applications of Linear Algebra in economic theory. To 
begin with, let us consider quantitative aspects of the exchange of goods between a certain number 
of economic agents. We here aim at a simplified abstract description of real economic processes. 


4.1 General considerations 

The quantitative model to be described in the following is due to the Russian 
economist Wassily Wassilyovich Leontief (1905-1999), cf. Leontief (1936) fl20l . 
for which, besides other important contributions, he was awarded the 1973 
Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel, 

Suppose given an economic system consisting of n G N interdependent economic agents ex¬ 
changing between them the goods they produce. For simplicity we want to assume that every one 
of these economic agents represents the production of a single good only. Presently we intend 
to monitor the flow of goods in this simple economic system during a specified reference period 
of time. The total numbers of the n goods leaving the production sector of this model constitute 
the OUTPUT quantities. The INPUT quantities to the production sector are twofold. On the 
one hand, there are exogenous INPUT quantities which we take to be given by m G N different 
kinds of external resources needed in differing proportions to produce the n goods. On the other 
hand, due to their mutual interdependence, some of the economic agents require goods made by 
tbeir neighbours to be able to produce their own goods; these constitute the endogenous INPUT 
quantities of the system. Likewise, the production sector’s total OUTPUT during the chosen ref¬ 
erence period of the n goods can be viewed to flow through one of two separate channels: (i) the 
exogenous channel linking the production sector to external consumers representing an open 
market, and (ii) the endogenous channel linking the economic agents to their neighbours (thus 
respresenting their interdependencies). It is supposed that momentum is injected into the economic 
system, triggering the flow of goods between the different actors, by the prospect of increasing 
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the value of the INPUT quantities, in line with the notion of the eeonomie value chain. 

Leontief’s model is based on the following three elementary 

Assumptions: 

1. For all goods involved the funetional relationship between INPUT and OUTPUT quantities 
be of a linear nature [cf. Eq. (12.51) 1. 

2. The proportions of “INPUT quantities to OUTPUT quantities” be constant over the refer- 
enee period of time considered; the flows of goods are thus considered to be stationary. 

3. Economic equilibrium obtains during the reference period of time: the numbers of goods 
then supplied equal the numbers of goods then demanded. 

The mathematical formulation of Leontief’s quantitative model employs the following 

Vector- and matrix-valued quantities: 

1. q — total output vector G components g* > 0 units 

2 . y — final demand vector G components > 0 units 

3. P — input-output matrix G components > 0 

4. (1 — P) — technology matrix G regular, hence, invertible 

5. (1 — P)~^ — total demand matrix G R"^^” 


6. V — resource vector G 


pmxl 


, components Vi > 0 units 


7. R — resource consumption matrix G 


\ components Rij > 0 , 


(dim: units) 
(dim: units) 
(dim: 1 ) 
(dim: 1 ) 
(dim: 1 ) 
(dim: units) 
(dim: 1 ) 


where 1 denotes the (n x n)-unit matrix [cf. Eq. (12.31) 1. Note that the components of all the 
vectors involved, as well as of the input-output matrix and of the resource consumption matrix, 
can assume non-negative values (!) only. 


4.2 Input-output matrix and resource consumption matrix 

We now turn to provide the definition of the two central matrix-valued quantities in Leontief’s 
model. We will also highlight their main characteristic features. 

4.2.1 Input-output matrix 

Suppose the reference period of time has ended for the economic system in question, i.e. the 
stationary flows of goods have stopped eventually. We now want to take stock of the numbers of 
goods that have been delivered by each of the n economic agents in the system. Say that during 
the reference period considered, agent 1 delivered of their good the number nn to themselves, the 
number 77-12 to agent 2, the number 7713 to agent 3, and so on, and, lastly, the number 77 i„ to agent n. 
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The number delivered by agent 1 to external eonsumers shall be denoted by r/i. Sinee in this model 
a good produeed cannot all of a sudden disappear again, and sinee by Assumption 3 above the 
number of goods supplied must be equal to the number of goods demanded, we find that for the 
total output of agent 1 it holds that qi := nu + ...+ ny + ... + uin + yi. Analogous relations 
hold for the total output q 2 , qs, ..., Qn of each of the remaining n — 1 agents. We thus obtain the 
intermediate result 


rill + ., 

• • + riij + . 

• • + riin + t/i > 0 

(4.1) 

na + .. 

• + rijj + .. 

.+nin + yi>0 

(4.2) 

Unl + . 

• • + ^nj + . 

• • T ri^jj A y^i > 0 . 

(4.3) 


This simple system of balance equations can be summarised in terms of a standard input-output 
table as follows: 


[Values in units] 

agent 1 ■■■ 

agent j ■ 

agent n 

external consumers 

S: total output 

agent 1 

rill 

rip 


yi 


agent i 

Tin 

Tlij 

• ^in 

yi 

Qi 

agent n 

TLnl 

Tlnj 

• ^nn 

yn 

Qn 


The first column of this table lists all the n different sources of flows of goods (or suppliers of 
goods), while its first row shows the n + 1 different sinks of flows of goods (or consumers of 
goods). The last column contains the total output of each of the n agents in the reference period 
of time. 

Next we compute for each of the n agents the respective values of the non-negative ratios 

INPUT (in units) of agent i for agent j (during reference period) 

OUTPUT (in units) of agent j (during reference period) 

or, employing a compact and economical index notation^ 


p_ 

^ij • — 


(4.5) 


'Note that the normalisation quantities in these ratios are given by the total output qj of the receiving agent j 
and not by the total output qi of the supplying agent i. In the latter case the P^ would represent percentages of the 
total output qi. 
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with i,j = 1,... ,n. These n x n = different ratios may be naturally viewed as the elements of 
a quadratie matrix P of format {n x n). In general, this matrix is given by 



/ nil 

nij 

nin \ 




nji-\-...-\-njj-\-...-\-njn-\-yj 

'Rnl + • • •+'Rri j + • • •+"^7171+yn 


T} 

nn 

riij 

nin 


Jr — 

nii + ...+nij+...+ni„+yi 

nji-\-...-\-njj-\-...-\-njn-\-yj 

Hnl + • • • + Rri j + • • • +'Rnn +^77 



i nnl 

n-nj 

Hnn 



\ n\\ + ...+n\j+...+nin+yi 

nji-\-...-\-njj-\-...-\-njn-\-yj 

'Rnl+---H“RnjH“---+'Rnn+yn / 



(4.6) 

and is referred to as Leontief’s input-output matrix of the stationary eeonomie system under 
investigation. 

For the very simple ease with just n = 3 produeing agents, the input-output matrix reduees to 


P 


nil 

ni2 

niz 

nii+ni2^nis^yi 

R21 

R21+R22+R23+?/2 

R22 

nzi +R32 +R33 +J/3 
R23 

nii+ni2+nis+yi 

nsi 

R21+R22+R23+?/2 

R32 

^31+^32+^33+J/3 

nss 

nii+ni2+ni3+yi 

n2i+n22+n23+y2 

n3i +n32 +n33 +y3 


It is important to realise that for an aetual eeonomie system the input-output matrix P ean be 
determined only once the reference period of time chosen has come to an end. 


The utility of Leontief’s stationary input-output matrix model is in its application for the purpose 
of forecasting. This is done on the basis of an extrapolation, namely by assuming that an input- 
output matrix Preference period Obtained from data taken during a specific reference period also is valid 
(to an acceptable degree of accuracy) during a subsequent period, i.e.. 



or, in component form. 


"^3 \ subsequent period 


^3 subsequent period 


I reference period 


^3 reference period 


In this way it becomes possible to compute for a given (idealised) economic system approximate 
numbers of INPUT quantities required during a near future production period from the known 
numbers of OUTPUT quantities of the most recent production period. Long-term empirical ex¬ 
perience has shown that this method generally leads to useful results to a reasonable approximation. 
All of these calculations are grounded on linear relationships describing the quantitative aspects of 
stationary flows of goods, as we will soon elucidate. 


4.2.2 Resource consumption matrix 

The second matrix-valued quantity central to Leontief’s stationary model is the resource con¬ 
sumption matrix R. This may be interpreted as providing a recipe for the amounts of the m dif¬ 
ferent kinds of external resources (the exogenous INPUT quantities) that are needed in the pro¬ 
duction of the n goods (the OUTPUT quantities). Its elements are defined as the ratios 

Rij ;= amounts (in units) required of resource i for the production of one unit of good j , 

(4.9) 
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with i = 1,... ,m und j = 1,..., n. The rows of matrix R thus contain information concerning the 
m resources, the columns information concerning the n goods. Note that in general the (m x n) 
resource consumption matrix R is not (!) a quadratic matrix and, therefore, in general not 
invertible. 


4.3 Stationary linear flows of goods 


4.3.1 Flows of goods: endogenous INPUT to total OUTPUT 


We now turn to a quantitative description of the stationary flows of goods that are associated 
with the total output q during a specific period of time considered. According to Leontief’s 
Assumption 1, there exists a linear functional relationship between the endogenous vector-valued 
INPUT quantity q — y and the vector-valued OUTPUT quantity q. This may be represented in 
terms of a matrix-valued relationship as 


q-y = Fq ^ q^-y^ = ^ P^-q^ , 

i=i 


(4.10) 


with i = 1,..., 77,, in which the input-output matrix P takes the role of mediating a mapping 
between either of these vector-valued quantities. According to Assumption 2, the elements of the 
input-output matrix P remain constant for the period of time considered, i.e. the corresponding 
flows of goods are assumed to be stationary. 

Relation (14.101) may also be motivated from an alternative perspective that takes the physical sci¬ 
ences as a guidline. Namely, the total numbers q of the n goods produced during the period of 
time considered which, by Assumption 3, are equal to the numbers supplied of the n goods satisfy 
a conservation law: “whatever has been produced of the n goods during the period of time con¬ 
sidered cannot get lost in this period.” In quantitative terms this simple relationship may be cast 
into the form 


total output final demand (exogenous) deliveries to production sector (endogenous) 

For computational purposes this central stationary flow of goods relation (14.101) may be rearranged 
as is convenient. In this context it is helpful to make use of the matrix identity q = Iq, where 1 
denotes the (n x n)-unit matrix [cf. Eq. (12.31) 1. 

Examples: 

(i) given/known: P, q 
Then it applies that 

n 

y = (1 - P)g yi = - Pij)(ij, (4.11) 

i=i 

with i = 1,... ,n; (1 — P) represents the invertible technology matrix of the economic 
system regarded. 
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(ii) given/known: P, y 
Then it holds that 


q = i^-P) ^ qi = - Pij) , (4.12) 

j=i 

with i = 1 ,..., n; (1 — P)“^ here denotes the total demand matrix, i.e., the inverse of the 
teehnology matrix. 

4.3.2 Flows of goods: exogenous INPUT to total OUTPUT 

Likewise, by Assumption 1, a linear funetional relationship is supposed to exist between the ex¬ 
ogenous veetor-valued INPUT quantity v and the veetor-valued OUTPUT quantity q. In matrix 
language this ean be expressed by 


u = Rq Vi = '^ RijQj , 


(4.13) 


with i = 1,..., m. By Assumption 2, the elements of the resource consumption matrix R remain 
constant during the period of time considered, i.e., the corresponding resource flows are supposed 

to be stationary. 


By combination of Eqs. (14.131) and (14.121) . it is possible to compute the numbers v of resources 
required (during the period of time considered) for the production of the n goods for given final 
demand y. It applies that 


V = Tlq = R(1 -P) ^y - Pjk) ^Vk , (4-14) 

i=l k=i 


with i = 1 ,..., m. 

GDC: For problems with n < 5, and known matrices P and R, Eqs. (14.1 II) . (14.121) and (14.141) can 
be immediately used to calculate the quantities q from given quantities y, or vice versa. 


4.4 Outlook 

Eeontief’s input-output matrix model may be extended in a straightforward fashion to include 
more advanced considerations of economic theory. Supposing a closed though not necessarily 
stationary economic system G comprising n interdependent economic agents producing n differ¬ 
ent goods, one may assign monetary values to the INPUT quantity v as well as to the OUTPUT 
quantities q and y of the system. Besides the numbers of goods produced and the associated 
flows of goods one may monitor with respect to G for a given period of time, one can in addition 
analyse in time and space the amount of money coupled to the different goods, and the corre¬ 
sponding flows of money. However, contrary to the number of goods, in general there does not 
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exist a conservation law for the amount of money with respeet to G. This may render the analysis 
of flows of money more diffieult, beeause, in the sense of an increase in value, money can either 
be generated inside G during the period of time considered or it can likewise be annihilated', it is 
not just limited to either flowing into respeetively flowing out of G. Central to considerations of 
this kind is a balance equation for the amount of money contained in G during a given period of 
time, which is an additive quantity. Such balance equations constitute familiar tools in Physics (cf. 
Herrmann (2003) [fTSl p 7ff]). Its structure in the present case is given bjH 


( rate of change in time 
of the amount of money 

in G [in CU/TU] 


flux of money \ f rate of generation of money 

into G [in CU/TU] in G [in CU/TU] 


Note that, with respect to G, both fluxes of money and rates of generation of money can in prin¬ 
ciple possess either sign, positive or negative. To deal with these quantitative issues properly, 
one requires the technical tools of the differential and integral calculus which we will discuss 
at an elementary level in Chs. |7]and [8] We make contact here with the interdisciplinary science 
of Econophysics (cf., e.g., Bouchaud and Potters (2003) [Q), a very interesting and challenging 
subject which, however, is beyond the scope of these lecture notes. 

Leontief’s input-output matrix model, and its possible extension as outlined here, provide the 
quantitative basis for considerations of economical ratios of the kind 

OUTPUT [in units] 

INPUT [in units] ’ 

as mentioned in the Introduction. In addition, dimensionless (scale-invariant) ratios of the form 

REVENUE [in CU] 

COSTS [in CU] ’ 

referred to as economic efficiency, can be computed for and compared between different economic 
systems and their underlying production sectors. In Ch.|7]we will briefly reconsider this issue. 


^Here the symbols CU and TU denote “currency units” and “time units,” respectively. 
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Chapter 5 

Linear programming 


On the backdrop of the economic principle, we discuss in this chapter a special class of quanti¬ 
tative problems that frequently arise in specific practical applications in Business and Manage¬ 
ment. Generally one distinguishes between two variants of the economic principle: either (i) to 
draw maximum utility from limited resources, or (ii) to reach a specific target with minimum ef¬ 
fort (costs). With regard to the ratio (OUTPUT)/(INPUT) put into focus in the Introduction, the 
issue is to find an optimal value for this ratio under given boundary conditions. This aim can be 
realised either (i) by increasing the (positive) value of the numerator for fixed (positive) value of 
the denominator, or (ii) by decreasing the (positive) value of the denominator for fixed (positive) 
value of the numerator. The class of quantitative problems to be looked at in some detail in this 
chapter typically relate to boundary conditions according to case (i). 


5.1 Exposition of a quantitative problem 


To be maximised is a (non-negative) real-valued quantity 2 ;, which depends in a linear functional 
fashion on a fixed number of n (non-negative) real-valued variables xi,..., We suppose that 
the n variables xi,..., Xn in turn are constrained by a fixed number m of algebraic conditions, 
which also are assumed to depend on xi,... ,Xn in a linear fashion. These m constraints, or 
restrictions, shall have the character of imposing upper limits on m different kinds of resources. 

Def.: Consider a matrix A G a vector b G two vectors c,x G and a constant 

d G M. A quantitative problem of the form 


max ■ X + d | Aa? < b, a; > 0 } , 


(5.1) 
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or, expressed in terms of a component notation, 


max z{xi, . . . , Xn) = CiXi + . . . + CnXn + d (5.2) 

auXi + ... + ainXn < bi (5.3) 

(^mlX\ . . . 3“ ClmnXn ^ (5.4) 

xi > 0 (5.5) 

Xn > 0, (5.6) 


is referred to as a standard maximum problem of linear programmin g with n real-valued vari¬ 
ables. The different quantities and relations appearing in this definition are called 


• z{xi ,..., Xn) — linear objective function, the dependent variable, 

• xi,... ,Xn — n independent variables, 

• Ax < b — m restrictions, 

• a; > 0 — n non-negativity constraints. 


Remark; In an analogous fashion one may also formulate a standard minimum problem of 
linear programming, which can be cast into the form 

min \^z = ■ X + d | Aa; > b, a; > 0 } . 


In this case, the components of the vector b need to be interpreted as lower limits on certain 
capacities. 


For given linear objective function z{xi,..., Xn), the set of points x 
the condition 

z{xi ,..., Xn) = C = constant G M , 


{xi,...,Xn)'^ satisfying 
(5.7) 


for fixed value of C, is referred to as an isoquant of 2:. Isoquants of linear objective functions 
of n = 2 independent variables constitute straight lines, of n = 3 independent variables Euclid¬ 
ian planes, of n = 4 independent variables Euclidian 3-spaces (or hyperplanes), and of n > 5 
independent variables Euclidian (n — l)-spaces (or hyperplanes). 


In the simplest cases of linear programmin g , the linear objective function 2 : depends on n = 2 
variables xi and X 2 only. An illustrative and efficient method of solving problems of this kind will 
be looked at in the following section. 
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5.2 Graphical method for solving problems with two indepen¬ 
dent variables 

The systematic graphical solution method of standard maximum problems of linear programming 
with n = 2 independent variables comprises the following steps: 

1 . Derivation of the linear objective function 

z{xi,X2) = CiXi + C 2 X 2 + d 

in dependence on the variables xi and X2. 

2. Identification in the xi, a: 2 -plane of the feasible region D of z which is determined by the 
m restrictions imposed on xi and X 2 . Specifically, D constitutes the domain of 2 ; (cf. Ch.|7]). 

3 . Plotting in the xi, X2-plane of the projection of the isoquant of the linear objective function 
which intersects the origin (0 = Xi = X2). When C2 7^ 0, this projection is described by the 
equation 

X 2 = -(Cl/C2)xi . 

4 . Erecting in the origin of the xi, X2-plane the direction of optimisation for z which is deter¬ 
mined by the constant ^-gradient 



5. Parallel displacement in the xi, X 2 -plane of the projection of the (0,0)-isoquant of 2 ; along 
the direction of optimisation (V^;)^ across the feasible region D out to a distance where the 
projected isoquant just about touches D. 

6 . Determination of the optimal solution (xio,X 2 o) as the point resp. set of points of inter¬ 
section between the displaced projection of the (0, 0)-isoquant of z and the far boundary 
ofD. 

7. Computation of the optimal value of the linear objective function zq = z{xio, X 20 ) from 
the optimal solution (xio, X 2 o)- 

8 . Specification of potential remaining resources by substitution of the optimal solution 
(xio, X 20 ) into the m restrictions. 

In general one finds that for a linear objective function z with n = 2 independent variables xi 
and X 2 , the feasible region D, when non-empty and bounded, constitutes an area in the xi,X 2 - 
plane with straight edges and a certain number of vertices. In these cases, the optimal values of 
the linear objective function 2 ; are always to be found either at the vertices or on the edges of the 
feasible region D. When D is an empty set, then there exists no solution to the corresponding 
linear programming problem. When D is unbounded, again there may not exist a solution to the 
linear programming problem, but this then depends on the specific circumstances that apply. 
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Remark: To solve a standard minimum problem of linear programming with n = 2 indepen¬ 
dent variables by means of the graphieal method, one needs to parallelly displaee in the Xi, X 2 - 
plane the projeetion of the (0, 0)-isoquant of 2 ; along the direetion of optimisation (until 
contact is made with the feasible region D for the first time. The optimal solution is then given by 
the point resp. set of points of intersection between the displaced projection of the (0, 0)-isoquant 
of z and the near boundary of D. 


5.3 Dantzig’s simplex algorithm 

The main disadvantage of the graphical solution method is its limitation to problems with 
only n = 2 independent variables. In actual practice, however, one is often concerned with 
linear programming problems that depend on more than two independent variables. To 
deal with these more complex problems in a systematic fashion, the US-American mathemati¬ 
cian George Bernard Dantzig (1914-2005) has devised during the 1940ies an efficient algorithm 
which can be programmed on a computer in a fairly straightforward fashion; cf. Dantzig 
(1949,1955) nil. 

In mathematics, simplex is an alternative name used to refer to a convex polyhedron, i.e., a body of 
finite (hyper-)volume in two or more dimensions bounded by linear (hyper-)surfaces which inter¬ 
sect in linear edges and vertices. In general the feasible regions of linear programming problems 
constitute such simplexes. Since the optimal solutions for the independent variables of linear 
programming problems, when they exist, are always to be found at a vertex or along an edge of 
simplex feasible regions, Dantzig developed his so-called simplex algorithm such that it system¬ 
atically scans the edges and vertices of a feasible region to identify the optimal solution (when it 
exists) in as few steps as possible. 

The starting point be a standard maximum problem of linear programming with n indepen¬ 
dent variables in the form of relations (I5.2I) - (I5.6I) . First, by introducing m non-negative slack 
variables si,..., Sm, one transforms the m linear restrictions (inequalities) into an equivalent set 
of m linear equations. In this way, potential differences between the left-hand and the right-hand 
sides of the m inequalities are represented by the slack variables. In combination with the defining 
equation of the linear objective function z, one thus is confronted with a system of 1 -I- m linear 
algebraic equations for thel + n + m variables z, xi, ..., si, ..., Sm, given by 


Maximum problem of linear programming in canonical form 

2 ; - CiXi - C 2 X 2 - ... - CnXn = d (5.8) 

aiiXi -f 012 X 2 -t-... -f- ainXn + Si = bi (5.9) 

O21X1 -f 022X2 -f . . . -[- Cl2nXn + S2 = &2 (5.10) 

T ®?n2^2 T . . . T (^mnXn T Sm • (5.11) 


As discussed previously in Ch. [3l a system of linear algebraic equations of format (1 + m) x 
(1 -f n -f m) is under-determined and so, at most, allows for multiple solutions. The general 
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(1 + n + m)-dimensional solution vector 

i^Ly ^l,Ly • • • y ^71,Ly ^l,Ly • • • y ^m,L') (5.12) 

thus contains n variables the values of which can be chosen arbitrarily. It is very important to be 
aware of this fact. It implies that, given the linear system is solvable in the first place, one has a 
choice amongst different solutions, and so one can pick the solution which proves optimal for the 
given problem at hand. Dantzig’s simplex algorithm constitues a tool for determining such an 
optimal solution in a systematic way. 

Let us begin by transferring the coefficients and right-hand sides (RHS) of the under-determined 
linear system introduced above into a particular kind of simplex tableau. 

Initial simplex tableau 


2 ; 

Xi 

X2 

Xn 

Si 

S2 

• • • 

RHS 

1 

-Cl 

-C2 

• • • C77, 

0 

0 

... 0 

d 

0 

Oil 

ai2 

• • • 

1 

0 

... 0 

hi 

0 

^21 

^22 

• • • 

0 

1 

... 0 

&2 

0 


^m2 

. . . Ojmn 

0 

0 

... 1 

hm 


In such a simplex tableau one distinguishes so-called basis variables from non-basis variables. 
Basis variables are those that contain in their respective columns in the number tableau a (1 -I- m)- 
component canonical unit vector [cf. Eq. (11.101) 1: in total the simplex tableau contains 1 -f m of 
these. Non-basis variables are the remaining ones that do not contain a canonical basis vector 
in their respective columns; there exist n of this kind. The complete basis can thus be perceived 
as spanning a (1 -f m)-dimensional Euclidian space Initially, always 2 : and the m slack 

variables si,..., constitute the basis variables, while the n independent variables xi,..., 
classify as non-basis variables [cf. the initial tableau (15.131) 1. The corresponding so-called (first) 
basis solution has the general appearance 

T T 

('^Si y ^l,Biy ■ ■ ■ y tCn,B\ y ) • • • ) ^m,Bi ) ((7,0,...,0,6i,..., &m) ; 

since, for simplicity, each of the n arbitrarily specifiable non-basis variables may be assigned the 
special value zero. In this respect basis solutions will always be special solutions (as opposed 
to general ones) of the under-determined system (I5.8I) - (I5.11I) — the maximum problem of linear 
programming in canonical form. 

Central aim of the simplex algorithm is to bring as many of the n independent vari¬ 
ables Xi,...,Xn as possible into the (1 -f m)-dimensional basis, at the expense of one of the 
m slack variables si,..., Sm, one at a time, in order to construct successively more favourable 
special vector-valued solutions to the optimisation problem at hand. Ultimately, the simplex algo- 
rit bm needs to be viewed as a special variant of GauBian elimination as discussed in Ch. [3l with a 
set of systematic instructions concerning allowable equivalence transformations of the underlying 
under-determined linear system (I5.8I) - (I5.11I) . resp. the initial simplex tableau (15.131) . This set of 
systematic algebraic simplex operations can be summarised as follows: 

Simplex operations 
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SI: Does the current simplex tableau show —Cj > 0 for all j G ,n}I If so, then the 

corresponding basis solution is optimal. END. Otherwise goto S2. 

S2: Choose a pivot column index j* e n} such that —Cj* := min{—Cj\j G 

n}} < 0. 

S3: Is there a row index i* G m} such that ai*j* > 0? If not, the objective function 2: is 

unbounded from above. END. Otherwise goto S4. 

S4: Choose a pivot row index such that 0 ^*^* > 0 and ;= mm{hi/ ai*j*\ai*j* > 0,i G 

m}}. Perform a pivot operation with the pivot element ai*j*. Goto SI. 

When the final simplex tableau has been arrived at, one again assigns the non-basis variables the 
value zero. The values of the final basis variables corresponding to the optimal solution of the 
given linear programmin g problem are then to be determined from the final simplex tableau by 
backward substitution, beginning at the bottom row. Note that slack variables with positive values 
belonging to the basis variables in the optimal solution provide immediate information on existing 
remaining capacities in the problem at hand. 



Chapter 6 

Elementary financial mathematics 


In this chapter we want to provide a brief introduetion into some basie eoneepts of financial math¬ 
ematics. As we will try to emphasise, many applieations of these eoneepts (that have immediate 
praetieal relevanee) are founded on only two simple and easily aeeessible mathematieal struetures: 
the so-called arithmetical and geometrical real-valued sequenees and their assoeiated finite series. 


6.1 Arithmetical and geometrical sequences and series 

6.1.1 Arithmetical sequence and series 

An arithmetical sequence of n G N real numbers a„ G M, 


(®n) 


nSN ) 


is defined by the property that the difference d between neighbouring elements in the sequence be 
constant, i.e., for n > 1 


On — ttn-i =: d = constant ^ 0 , 


( 6 . 1 ) 


with a„, a„_i, d G M. Given this reeursive formation rule, one may infer the explicit representa¬ 
tion of an arithmetical sequence as 


On = ai + {n — l)d with n G M . (6.2) 

Note that any arithmetical sequence is uniquely determined by the two free parameters Oi and d, 
the starting value of the sequenee and the eonstant differenee between neighbours in the sequence, 
respeetively. Equation (16.21) shows that the elements in a non-trivial arithmetical sequence 
exhibit either linear growth or linear decay with n. 

When one calculates for an arithmetical sequence of n -I-1 real numbers the arithmetical mean 
of the immediate neighbours of any particular element a„ (with n > 2), one finds that 

- {un-i -f Un+i) = 2 ~ + Ui + nd) = ai + {n — l)d = a„ . (6.3) 
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Summation of the first n elements of an arbitrary arithmetical sequence of real numbers leads to 

a finite arithmetical series, 


d 


S'n := Oi + 02 + ... + a„ = ^ Ofc = ^ [oi + (fc - l)d] = nai + - {n - l)n . 


(6.4) 


k=l 


k=l 


In the last algebraie step use was made of the GauBian identitjQ(cf., e.g., Boseh (2003) |I3 p 21]) 


It — j. 

'^k = -{n-l)n. 

k=l 


(6.5) 


6.1.2 Geometrical sequence and series 


A geometrical sequence of n e N real numbers a„ e R, 

(On)n£N 5 


is defined by the property that the quotient q between neighbouring elements in the sequenee be 
constant, i.e., for n > 1 


—— ='. q = constant ^ 0 , 

O-n—1 


( 6 . 6 ) 


with anyttn-i G R and q G R\{0,1}. Given this recursive formation rule, one may infer the 

explicit representation of a geometrical sequence as 


On = oig” ^ with n G M. 


(6.7) 


Note that any geometrical sequence is uniquely determined by the two free parameters oi and q, 
the starting value of the sequence and the constant quotient between neighbours in the sequence, 
respectively. Equation (16.71) shows that the elements in a non-trivial geometrical sequence 
exhibit either exponential growth or exponential decay with n (cf. Sec. 17.1.41) . 

When one calculates for a geometrical sequence of n + 1 real numbers the geometrical mean of 
the immediate neighbours of any particular element a„ (with n > 2), one finds that 

^Un-i ■ Un+i = s/aiq^-'^ ■ aiq^ = . ( 6 . 8 ) 


Summation of the first n elements of an arbitrary geometrical sequence of real numbers leads to 

a finite geometrical series. 


n—1 


Sn ■= tti + 02 + ■ ■ ■ + On = Ufc = [oiq^ = Oi q^ = ai ——^ . 

k=l k=l k=0 ^ 

In the last algebraic step use was made of the identity (cf., e.g., Bosch (2003) |I3 p 27]) 


(6.9) 



( 6 . 10 ) 


n 

'Analogously, the modified GauBian identity — 1) = applies. 

k^l 
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6.2 Interest and compound interest 


Let us consider a first rather simple interest model. Suppose given an initial capital of positive 
value Kq > 0 CU paid into a bank account at some initial instant, and a time interval consisting 
of n G N periods of equal lengths. At the end of each period, the money in this bank account 
shall earn a service fee corresponding to an interest rate of p > 0 percent. Introducing the 
dimensionless interest factoiH 

+ ( 6 . 11 ) 

one finds that by the end of the first interest period a total capital of value (in CU) 


K, 


Ko + Ko 


P 

100 


= Ko I + 


P \ ^ 
1007 


Koq 


will have accumulated. When the entire time interval of n interest periods has ended, a final 
capital worth of (in CU) 


recursively: Kn = Kn-iq , n G N , 


( 6 . 12 ) 


will have accumulated, where Kn-i denotes the capital (in CU) accumulated by the end of n — 1 
interest periods. This recursive representation of the growth of the initial capital Kq due to a total 
of n interest payments and the effect of compound interest makes explicit the direct link with the 
mathematical structure of a geometrical sequence of real numbers (16.61) . 

It is a straightforward exercise to show that in this simple interest model the final capital Kn is 
related to the initial capital Kq by 


explicitly: = A'oq'" , n G N . 


(6.13) 


Note that this equation links the four non-negative quantities Kn, Kq, q and n to one another. 
Hence, knowing the values of three of these quantities, one may solve Eq. (16.131) to obtain the 
value of the fourth. For example, solving Eq. (16.131) for Kq yields 

Kq = ^=-.Bq. (6.14) 

qn 

In this particular variant, Kq is referred to as the present value Bq of the final capital Kn', this is 
obtained from Kn by an ?7,-fold division with the interest factor q. 

Further possibilities of re-arranging Eq. (16.131) are: 


(i) Solving for the interest factor q: 



(ii) Solving for the contract period n: 


\n{Kn/KQ) 

\n{q) 


(6.15) 


(6.16) 


^Inverting this defining relation for q leads to p = 100 • (q — 1). 
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From now on, n G N shall denote the number of full years that have passed in a speeifie interest 
model. 

Now we turn to diseuss a seeond, more refined interest model. Let us suppose that an initial 
capital > 0 CU earns interest during one full year m G N times at the mth part of a nominal 
annual interest rate pnom > 0. At the end of the first out of m periods of equal length 1/m, the 
initial eapital Kq will thus have inereased to an amount 


Kl/m — Kq + Kq 


Pn 


= Kn 1 + 


Pi, 




m- 100 V ' m- 100/ 

By the end of the kth (k < m) out of m periods the account balance will have beeome 


the interest faetor (1 + 


Pnom \ 

) will then have been applied k times to Kq. At the end of the full 


m ■ 100 / 

year, Kq in this interest model will have inereased to 

= K^im = Ko(l + , m G N . 

This relation defines an effective interest factor 

/, , Pnom 

QeS :— ( 1 /-77^ ) , 

V m■100/ 

with associated effective annual interest rate 


(6.17) 


Peff = 100 ■ 



Pnom 

m ■ 100 



m G N , 


(6.18) 


obtained from re-arranging = 1 + 

When, ultimately, n G N full years will have passed in the second interest model, the initial 
capital Kq will have been transformed into a final capital of value 


Kn 


Ko 



Pnom 

m ■ 100 


n-m 


Koq^Q , n, m G N . 


(6.19) 


The present value Bq of is thus given by 


Bo = ^ = A-„ 

7eff 


( 6 . 20 ) 


Finally, as a third interest model relevant to applications in Finance, we turn to consider the con¬ 
cept of installment savings. For simplicity, let us restrict our discussion to the case when n G N 
equal installments of constant value E > 0 CU are paid into an account that earns p > 0 percent 
annual interest (i.e., g > 1) at the beginning of each of n full years. The initial account balance 
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be i^o = 0 CU. At the end of a first full year in this interest model, the aceount balance will have 
increased to 

K^ = E + E-^ = e(i + -^) = Eq . 

100 V 100/ ^ 

At the end of two full years one finds, substituting for Ki, 


K2 = {Ki + E)q = {Eq + E)q = E{q^ + q) = Eq{q + 1) . 

At the end of n full years we have, recursively substituting for Kn-i, Kn-2, etc., 

n —1 

+ E)q = --- = E{q^ + . . . + + q) = Eq{q^-^ + . . . + g + l) = ^ . 

k=0 

Using the identity (16.101) . since presently g > 1, the account balance at the end of n full years can 
be reduced to the expression 


gn _ 1 

Kn = Eq^ - 

g G M>i , 

n e N . 

g- 1 




( 6 . 21 ) 


The present value Bq associated with is obtained by n-fold 
factor q\ _ 

^ E{q^ - 1 ) 

0- qu gn-l('g_X) 


division of Kn with the interest 


( 6 . 22 ) 


This gives the value of an initial capital Bq which will grow to the same final value Kn after n 
annual interest periods with constant interest factor g > 1. 


Lastly, re-arranging Eq. (16.211) to solve for the contract period n yields. 


\n[lE{q-l){Kn/Eq)] 

ln(g) 


(6.23) 


6.3 Redemption payments in constant annuities 

The starting point of the next discussion be a mortgage loan of amount > 0 CU that an 
economic agent borrowed from a bank at the obligation of annual service payments of p > 0 
percent (i.e., g > 1) on the remaining debt. We suppose that the contract between the agent and 
the bank fixes the following conditions: 

(i) the first redemption payment Ti amount to f > 0 percent of the mortgage Rq, 

(ii) the remaining debt shall be paid back to the bank in constant annuities of value A > 0 CU 
at the end of each full year that has passed. 

The annuity A is defined as the sum of the variable nth interest payment Zn > 0 CU and the 
variable nth redemption payment T„ > 0 CU. In the present model we impose on the annuity 
the condition that it be constant across full years, 

A = Zn + Tn = constant. 


(6.24) 
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For n = 1, for example, we thus obtain 


A = Zi + Ti = Ro ■ — + Ro ■ — = Ro 

“ 100 100 


p + t 
100 


— Rn 


(g- 1) + 


100 


= constant. (6.25) 


For the first full year of a running mortgage contract, the interest payment, the redemption payment, 
and, following the payment of a first annuity, the remaining debt take the values 

= -Ro ■ ^ = Ro{q - 1) 

Ti = A-Zi 

substitute for 2 

Ri = Rq + Zi — A 

By the end of a second full year, these become 

Z2 = Ri{q- 1 ) 

T2 = A — Z2 

R2 = Ri + Z2 — A 


Rq + Rq ' — A — Rqq — A . 


substitute for Z 2 


substitute for Ri 


Riq — A 


Roq^ ~ ^{q + 1 ) • 


At this stage, it has become clear according to which patterns the different quantities involved in 
the redemption payment model need to be formed. The interest payment for the nth full year in a 
mortgage contract of constant anuities amounts to (recursively) 

Zn = Rn-i{q - I) , neN, (6.26) 


where Rn-i denotes the remaining debt at the end of the previous full year. The redemption 
payment for full year n is then given by (recursively) 


Tn = A - Zn , n e N . 

The remaining debt at the end of the nth full year then is (in CU) 

recursively: Rn = Rn-i + Zn — A = Rn-iq — A , n G N . 

By successive backward substitution for Rn- 2 , etc., Rn can be re-expressed as 


n—1 

Rn = Roq^ ~ A(q^ ^ + ... + g + 1) = i?og” ~ A g^ . 

k=0 

Now employing the identity (16.101) . we finally obtain (since g > 1) 


explicitly: Rn 


g" - 1 

Roq^ - A ^ -- , n e N . 

g - 1 


(6.27) 


(6.28) 


(6.29) 


All the formulae we have now derived for computing the values of the quantities {n, T„, Rn} 

form the basis of a formal redemption payment plan, given by 
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n 

[CU] 

Tn [CU] 

Rn [CU] 

0 

- 

- 

Ro 

1 

Zi 

Ti 

Ri 

2 

Z2 

T2 

R2 


a standard scheme that banks must make available to their mortgage eustomers for the purpose of 
finaneial orientation. 

Remark: For known values of the free parameters > 0 CU, g > 1 and A > 0 CU, the simple 
reeursive formulae (16.261) . (16.271) and (16.281) ean be used to implement a redemption payment plan 
in a modern spreadsheet programme sueh as EXCEL or OpenOffiee. 

We emphasise the following observation eoneeming Eq. (16.291) : sinee the eonstant annuity A eon- 
tains implieitly a faetor (g — 1) [of. Eq. (16.251) 1. the two oompeting terms in this relation eaeh grow 
exponentially with n. Eor the redemption payments to eventually terminate, it is thus essential to 
fix the free parameter t (for known p > 0 <(=)> g > 1) in sueh a way that the seoond term on the 
right-hand side of Eq. (16.291) is given the possibility to oateh up with the first as n progresses (the 
latter of whioh has a head start of i?o > 0 CU at n = 0). The necessary condition following from 

I 

the requirement that < Rn-i is thus f > 0. 

Equation (16.291) links the five non-negative quantities Rn, Ro, g, n and A to one another. Given 
one knows the values of four of these, one can solve for the fifth. Eor example: 


(i) Calculation of the contract period n of a mortgage eontract, knowing the mortgage Rq, the 
interest faetor g and the annuity A. Solving the eondition = 0 imposed on R^ for n 
yields (after a few algebraie steps) 


n = 


1° (1 + f) 

ln(?) 




the eontract period is thus independent of the value of the mortgage loan, Rq. 


(6.30) 


(ii) Evaluation of the annuity A, knowing the eontraet period n, the mortgage loan Rq, and the 
interest faetor g. Solving the eondition Rn = 0 imposed on Rn for A immediately yields 


A 


g” -1 


Ro ■ 


(6.31) 


Now equating the two expressions (16.311) and (16.251) for the annuity A, one finds in addition 
that 

t g — 1 

Too “ g’^ - 1 ■ 


(6.32) 
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6.4 Pension calculations 


Quantitative models for pension calculations assume given an initial capital Kq > 0 CU that 
was paid into a bank account at a particular moment in time. The issue is to monitor the subse¬ 
quent evolution in discrete time n of the account balance Kn (in CU), which is subjected to two 
competing influences: on the one-hand side, the bank account earns interest at an annual interest 
rate of p > 0 percent (i.e., q > 1), on the other, it is supposed that throughout one full year a total 
of m G N pension payments of the constant amount a are made from this bank account, always 
at the beginning of each of m intervals of equal duration per year. 

Let us begin by evaluating the amount of interest earned per year by the bank account. An impor¬ 
tant point in this respect is the fact that throughout one full year there is a total of m deductions of 
value a from the bank account, i.e., in general the account balance does not stay constant through¬ 
out that year but rather decreases in discrete steps. For this reason, the account is credited by the 
bank with interest only at the mth part of p > 0 percent for each interval (out of the total of m) 
that has passed, with no compound interest effect. Hence, at the end of the first out of m intervals 
per year the bank account has earned interest worth of (in CU) 

Zi/„ = (A'„ - a) . = (A„ - a) LUf . 

m ■lUU m 

The interest earned for the kth interval (out of m; k < m) is then given by 

Zk/m = {Ko - ka) . 

m 

Summation over the contributions of each of the m intervals to the interest earned then yields for 
the entire interest earned during the first full year (in CU) 


Zi = ^ Zk/m = - ka) 


(g- 1) (g- 1) 


k=l 


k=l 


m 


m 


uiRq — a k 


k=l 


By means of substitution from the identity (16.5L this result can be recast into the equivalent form 


Zi = 


Kq- -{m+ l)a 


(9 - 1 ) • 


(6.33) 


Note that this quantity decreases linearly with the number of deductions m made per year resp. with 
the pension payment amount a. 

One now finds that the account balance at the end of the first full year that has passed is given by 


Eq. 

Ki = Kq — ma + Zi K^q 


m + -{m+ l)(g - 1) 


a . 


At the end of a second full year of the pension payment contract the interest earned is 


Zo = 


Ki - -{m + l)a 


(9 - 1 ) , 
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while the aeeount balanee amounts to 

substitute for Ki and Z2 

K 2 = Ki- ma + Z 2 


Koq^ 


m + - (m + l)(g — 1) 


a{q + 1) . 


At this stage, eertain fairly simple patterns for the interest earned during full year n, and the 
account balance after n full years, reveal themselves. For Zn we have 


Zr, — 


Rn-1 - l)a 


(7 - 1) , 


(6.34) 


and for K„ one obtains 


Rn R-n—l 


ma 


substitute for Kn—i and Zn 

Zn ^ 


m + - (m + l)(g - 1) 


n—1 

k=0 


The latter result ean be re-expressed upon substitution from the identity (16.101) . Thus, can 
finally be given by 


explieitly: Kn = Kgq^ — 


m + -(m + l)(q - 1) 


g” - 1 

a -, n,m E N 

g - 1 


(6.35) 


In a fashion praetieally identieal to our discussion of the redemption payment model in Sec. 16.31 
the two competing terms on the right-hand side of Eq. (16.351) likewise exhibit exponential growth 
with the number n of full years passed. Specifically, it depends on the values of the parameters 
Ko > 0 CU, g > 1, a > 0 CU, as well as m > 1, whether the second term eventually manages to 
catch up with the first as n progresses (the latter of which, in this model, is given a head start of 
value Ko > 0 CU at n = 0). 

We remark that Eq. (16.351) . again, may be algebraically re-arranged at one’s convenience (as long 
as division by zero is avoided). Eor example: 


(i) The duration n (in full years) of a particular pension contract is obtained from solving the 
condition Kn = 0 accordingly. Given that [.. .]a — Ko{q — 1) > 0, one thus find^ 


n = 


In 


( _L>_"i 

\l...]a-Ko{q-l) ) 
ln(g) 


(6.36) 


(ii) The present value Bo of a pension scheme results from the following consideration: for 
fixed interest factor g > 1, which initial capital iCo > 0 CU must be paid into a bank account 
such that for a duration of n full years one can receive payments of constant amount a at the 
beginning of each of m intervals (of equal length) per year? The value of Bo = Ko is again 
obtained from imposing on Eq. (16.351) the condition Kn = 0 and solving for Ko. This yields 


Bn = Kn = 


m + -{m + l){q — 1) 


g^(g- 1) 


(6.37) 


^To avoid notational overload, the brackets [...] here represent the term \Tn + ^ {m + l){q — 1)]. 
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(iii) The idea of so-ealled everlasting pension payments of amount Oever > 0 CU is based on the 
strategy to eonsume only the annual interest earned by an initial eapital Kq > 0 CU residing 
in a bank aeeount with interest faetor g > 1. Imposing now on Eq. (16.351) the eondition 

I 

Kn = Kq to hold for all values of n, and then solving for a, yields the result 




m + i (m + l)(g + 1) 


Kn 


(6.38) 


Note that, naturally, Oever is direetly proportional to the initial eapital Kq ! 


6.5 Linear and declining-balance depreciation methods 

Attempts at the quantitative deseription of the proeess of deelining material value of industrial 
goods, properties or other assets, of initial value Kq > 0 CU, are referred to as depreciation. In¬ 
ternational tax laws generally provide investors with a choice between two particular mathematical 
methods of calculating depreciation. We will discuss these options in turn. 


6.5.1 Linear depreciation method 

When the initial value Kq > 0 CU is supposed to decline to 0 CU in the space of N full years by 
equal annual amounts, the remaining value Rn (in CU) at the end of n full years is described by 


Rn Kq Tl 


\N 


n = 1,..., N . 


(6.39) 


Note that for the difference of remaining values for years adjacent one obtains Rn — Rn-i = 

=: d < 0. The underlying mathematical structure of the straight line depreciation 

method is thus an arithmetical sequence of real numbers, with constant negative difference d 
between neighbouring elements (cf. Sec. 16.1.11) . 



6.5.2 Declining-balance depreciation method 


The foundation of the second depreciation method to be described here, for an industrial good of 
initial value Kq > 0 CU, is the idea that per year the value declines by a certain percentage rate 
p > 0 of the value of the good during the previous year. Introducing a dimensionless depreciation 
factor by 


q:=l 


P 

100 


< 1 


(6.40) 


the remaining value Rn (in CU) after n full years amounts to 


recursively: Rn = Rn-iQ , Ro = Kq , n eN . 


(6.41) 


The underlying mathematical structure of the declining balance depreciation method is thus a 
geometrical sequence of real numbers, with constant ratio 0 < g < 1 between neighbouring 
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elements (ef. See. l6.1.2l) . With inereasing n the values of these elements beeome ever smaller. By 
means of sueeessive baekward substitution expression (16.411) can be transformed to 


explicitly: = Kqq'^ , 0<g<l, nGN. 


(6.42) 


From Eq. (16.421) . one may derive results concerning the following questions of a quantitative na¬ 
ture: 


(i) Suppose given a depreciation factor q and a projected remaining value Rn for some industrial 
good. After which depreciation period n will this value be attained? One finds 


InjRjKo) 

\n{q) 


(6.43) 


(ii) Knowing a projected depreciation period n and corresponding remaining value Rn, at which 
percentage rate p > 0 must the depreciation method be operated? This yields 


q = 






(6.44) 


6.6 Summarising formula 


To conclude this chapter, let us summarise the results on elementary financial mathematics that 
we derived along the way. Remarkably, these can be condensed in a single formula which contains 
the different concepts discussed as special cases. This formula, in which n represents the number 
of full years that have passed, is given by (cf. Zeh-Marschke (2010) ll26l l: 


Kn 


Koq^ + R 


qn _ I 

g- 1 ’ 


gGM>o\{l}, neN. 


(6.45) 


The different special cases contained therein are: 

(i) Compound interest for an initial capital Kq > 0 CU: with R = 0 and g > 1, Eq. (16.451) 
reduces to Eq. (16.131) . 

(ii) Installment savings with constant installments E > 0 CU: with Kq = 0 CU, g > 1 and 
R = Eq, Eq. (16.451) reduces to Eq. (16.191) . 

(iii) Redemption payments in constant annuities: with Kq = —Rq < 0 CU, g > 1 and 

R = A> 0 CU, Eq. (16.451) reduces to the negative (!) of Eq. (16.291) . In this dual formulation, 
remaining debt Kn = —Rn is (meaningfully) expressed as a negative account balance. 


(iv) Pension payments of constant amount a 

— m H—(m + l)(g —1) 


> 0 CU: with g > 1 and R = 

a, Eq. (16.451) transforms to Eq. (16.351) . 


(v) Declining balance depreciation of an asset of initial value Kq > 0 CU: with R = 0 and 
0 < g < 1, Eq. (16.451) converts to Eq. (16.421) for the remaining value Kn = Rn- 
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Chapter 7 

Differential calculus of real-valued 
functions of one real variable 


In Chs. [U to [5] of these leeture notes, we eonfined our eonsiderations to funetional relationships 
between INPUT quantities and OUTPUT quantities of a linear nature. In this ehapter now, we 
turn to diseuss eharaeteristic properties of truly non-linear functional relationships between one 

INPUT quantity and one OUTPUT quantity. 


7.1 Real-valued functions of one real variable 


Let us begin by introdueing the eoneeM of a real-valued function of one real variable. This 
constitutes a special kind of a mapping that needs to satisfy the following simple but very strict 
rule: 


a mapping / that assigns to every element x from a subset D of the real numbers M (i.e., 
D CW) one and only one element y from a second subset W of the real numbers M (i.e., W C M). 


Def.: A unique mapping / of a subset D C R of the real numbers onto a subset IL C R of the 
real numbers. 


f-.D^W, xv^y = f{x) 


(7.1) 


is referred to as a real-valued function of one real variable. 


We now fix some terminology concerning the concept of a real-valued function of one real variable: 


• D: domain of /, 


• W: target space of /, 


• X e D: independent variable of /, also referred to as the argument of /, 


• y eW: dependent variable of /, 

*Cf. our introduction in Ch.|2]of matrices as a particular class of mathematical objects. 
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• /(x): mapping prescription, 

• graph of /: the set of pairs of values G = {{x, f{x))\x e D} C 

For later analysis of the mathematieal properties of real-valued funetions of one real variable, we 
need to address a few more teehnieal issues. 

Def.: Given a mapping / that is one-to-one and onto, with domain D{f) C R and target spaee 
if) ^ not only is every x G Dif) assigned to one and only one y G Wif), but also every 
y G Wif) is assigned to one and only one x G Dif). In this ease, there exists an assoeiated 
mapping with Dif~^) = Wif) and Wif~^) = Dif), whieh is referred to as the inverse 
function of /. 

Def.: A real-valued funetion / of one real variable x is continuous at some value x G Dif) when 
for Ax G R>o the eondition 

lim fix — Ax) = lim fix -f Ax) = fix) (7.2) 

obtains, i.e., when at x the left and right limits of the funetion / eoineide and are equal to the value 
fix). A real-valued funetion / as sueh is continuous when / is eontinuous/or all x G Dif). 

Def.: When a real-valued funetion / of one real variable x satisfies the eondition 

/(a) </( 6 ) for all a, 6 G £>(/) with a < 6 , (7.3) 

then / is ealled strictly monotonously increasing. When, however, / satisfies the eondition 

/(a) > fib) for all a,b G Dif) with a < b , (7.4) 

then / is ealled strictly monotonously decreasing. 

Note, in partieular, that real-valued funetions of one real variable that are strietly monotonous and 
eontinous are always one-to-one and onto and, therefore, are invertible. 

In the following, we briefly review five elementary elasses of real-valued funetions of one real vari¬ 
able that find frequent applieation in the modelling of quantitative problems in economic theory. 

7.1.1 Polynomials of degree n 

Polynomials of degree n are real-valued funetions of one real variable of the form 


y = fix) = OnX"^ -f ttn-ix"^ ^ -f ... -f ttiX^ -f ... -f a 2 X^ -\- aix + ao 
with Oj G R, f = 1,..., n, n G N, a„ 7 ^ 0 . 


(7.5) 


Their domain eomprises the entire set of real numbers, i.e., Dif) = R. The extent of their target 
spaee depends speeifieally on the values of the real eonstant coefficients a* G R. Functions in this 
class possess a maximum of n real roots. 
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7.1.2 Rational functions 


Rational functions are constructed by forming the ratio of two polynomials of degrees m resp. n. 


i.e., 


y = fix) 


Pmjx) 

qnix) 


ttmX'^ + . . . + aiX + Oo 
hnX'^ + .. . + hiX + hQ 


with Oj, bj G M, 


1,..., m, j = 1,..., n, m,n eN, am,bn^0 ■ 


(7.6) 


Their domain is given by D{f) = R\{x|g„(a:) = 0}. When for the degrees m and n of the 
polynomials Pmix) and g„(x) we have 


(i) m < n, then / is referred to as a proper rational function, or 


(ii) m> n, then / is referred to as an improper rational function. 

In the latter case, application of polynomial division leads to a separation of / into a purely poly¬ 
nomial part and a proper rational part. The roots of / always eorrespond to those roots of the 
numerator polynomial Pmix) for whieh simultaneously qnix) ^ 0 applies. The roots of the de¬ 
nominator polynomial qnix) eonstitute poles of /. Proper rational funetions always tend for very 
small (i.e., x —)■ —cxd) and for very large (i.e., x —)■ -|-oo) values of their argument to zero. 


7.1.3 Power-law functions 


Power-law functions exhibit the speeific structure given by 

y = fix) = a:" with a G M . 


(7.7) 


We here confine ourselves to oases with domains D{f) = M>o, suoh that for the ooresponding 
target spaces we have Wif) = M>o. Under these conditions, power-law functions are striotly 
monotonously inoreasing when a > 0, and strictly monotonously decreasing when a < 0. Henoe, 
they are inverted hy y = ^fx = There do not exist any roots under the oonditions stated here. 


7.1.4 Exponential functions 


Exponential functions have the general form 

y = fix) = with a G M>o\{l} . 


(7.8) 


Their domain is D{f) = M, while their target space is Wif) = M>o. They exhibit strict 
monotonous increase for a > 1, and striet monotonous deerease for 0 < a < 1. Henee, they 
too are invertible. Their j/-intereept is generally loeated aty = 1. For a > 1, exponential funetions 
are also known as growth functions. 


Speeial ease: When the constant (!) base number is ehosen to be a = e, where e denotes the 
irrational Euler’s number (according to the Swiss mathematician Leonhard Euler, 1707-1783) 
defined by the infinite series 


e 


CXD 


E 


1 


1 

0 ! 


+ 


1 

1 ! 


+ 


1 


+ 


1 


+ 


• • ? 
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one obtains the natural exponential function 

y = f\x) = =: exp(a;) . 

In analogy to the definition of e, the relation 


(7.9) 


= exp(a;) = 


0 19 *^ 

00^ 00^ 0C'~' 


oo h 
X 

k=0 


applies. 

7.1.5 Logarithmic functions 

Logarithmic functions, denoted by 


y = fix) = log„(a;) with a G M>o\{l} , 


(7.10) 


are defined as inverse functions of the strictly monotonous exponential functions y = f{x) = 

— and vice versa. Correspondingly, D{f) = R>o and Wif) = M apply. Strictly monotonously 
increasing behaviour is given when a > 1, strictly monotonously decreasing behaviour when 
0 < a < 1. In general, the x-intercept is located at a; = 1. 

Special case: The natural logarithmic function (lat.: logarithmus naturalis) obtains when the 
constant basis number is set to a = e. This yields 


y = fix) = loge(a;) ;= ln(a;) . 


(7.11) 


7.1.6 Concatenations of real-valued functions 

Real-valued functions from all five categories considered in the previous sections may be combined 
arbitrarily (respecting relevant computational rules), either via the four fundamental arithmetical 
operations, or via concatenations. 

Theorem; Let real-valued functions / and g be continuous on domains D{f) resp. D{g). Then 
the combined real-valued functions 

1. sum/difference f ± g, where (/ ± g)ix) := fix) ± gix) with D{f) n D{g), 

2. product / ■ g, where (/ ■ g)ix) ;= fix)gix) with D{f ) n D{g), 

3. quotient —, where (ix) ■= with g{x) 7 ^ 0 and D{f) fl D{g)\{x\gix) = 0}, 

9 \gj gix) 

4. concatenation fog, where (/ o g)ix) := fig{x)) mit {x G D{g)\g{x) G D{f)}, 

are also continuous on the respective domains. 
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7.2 Derivation of differentiable real-valued functions 
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The central theme of this chapter is the mathematical description of the local variability of con¬ 
tinuous real-valued function of one real variable, / iDCM—)-lVCM. To this end, let 
us consider the effect on / of a small change of its argument x. Supposing we affect a change 
a; —)■ a; + Aa;, with Aa; G M, what are the resultant consequences for /? We immediately find that 
y ^ y -\- Ay = f{x + Ax), with Ay G M, obtains. Hence, a prescribed change of the argument x 
by a (small) value Ax triggers in / a change by the amount Ay = f{x + Ax) — f{x). It is of 
general quantitative interest to compare the sizes of these two changes. This is accomplished by 
forming the respective difference quotient 

^ f{x + Aa:) - f{x) 

Ax Ax 

It is then natural, for given /, to investigate the limit behaviour of this difference quotient as the 
change Aa: of the argument of / is made successively smaller. 

Def.: A continuous real-valued function / of one real variable is called differentiable at 
X E D{f) , when for arbitrary Aa: G M the limit 


fix) 


1 . 

lim —— 

Aa;->-0 Ax 


lim 

Ax—>0 


/(x + Ax) - /(x) 
Ax 


(7.12) 


exists and is unique. When / is differentiable/or all x G D{f), then / as such is referred to as 
being differentiable. 

The existence of this limit in a point (x, /(x)) for a real-valued function / requires that the latter 
exhibits neither “jumps” nor “kinks,” i.e., that at (x,/(x)) the function is sufficiently “smooth.” 
The quantity /'(x) is referred to as the first derivative of the (differentiable) function / at posi¬ 
tion X. It provides a quantitative measure for the local rate of change of the function / in the 
point (x, fix)). In general one interprets the first derivative fix) as follows: an increase of the 
argument x of a differentiable real-valued function / by 1 (one) unit leads to a change in the value 
of / by approximately fix) ■ 1 units. 

Alternative notation for the first derivative of /: 


fix) 


dfjx) 

dx 


The differential calculus was developed in parallel with the integral calculus (see Ch.|7]) during the 
second half of the 17*^ Century, independent of one another by the English physicist, mathemat- 
iccian, astronomer and philosopher Sir Isaac Newton (1643-1727) and the German philosopher, 
mathematician and physicist Gottfried Wilhelm Leibniz (1646-1716), 

Via the first derivative of a differentiable function / at an argument xq G D(f), i.e., /'(xq), one 
defines the so-called linearisation of / in a neighbourhood of xq. The equation describing the 
associated tangent to / in the point (xo, fixo)) is given by 


y = fixo) + f ixo)ix - Xo) . 


(7.13) 
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GDC: Local values f'{xo) of first derivatives ean be eomputed for given funetion / in mode CALC 
using the interaetive routine dy/dx. 

The following rules of differentiation apply for the five families of elementary real-valued funetions 
diseussed in See. l7.li as well as eoneatenations thereof: 


Rules of differentiation 

1. (c)' = 0 for c = eonstant G M 

2 . (x)' = 1 

3. (x"')' = nx‘^~^ for n G N 

4. [x^y = ax°‘~^ for a G M and x G M>o 

5. {a^y = ln(a)a^ for a G M>o\{l} 

6. {e°'^y = ae°'^ for a G M 


7- (loga(a;))' = 


X ln(a^ 


for a G M>o\{l} and x G M>o 


8. (ln(a;))' = — forx G M>o 

X 

For differentiable real-valued functions / and g it holds that: 
1. {cf{x)y = cf'{x) for c = constant G M 
2- {f{x) ± g{x)y = f{x) ± g’{x) 

3. if{x)g{x)y = f{x)g{x) + f{x)g'{x) 

7(x) \ ' f{x)g{x) - f{x)g'{x) 


n{x)J {9{x)y 

5. {{fog){x)y=f{g)\^^^^^ygyx) 

6. (ln(/(a:)))' = for /(x) > 0 


V. {f-\x)y = 


fix) 
1 


f'iy) 


y=f 


(constants) 
(linear function) 
(natural power-law functions) 
(general power-law functions) 
(exponential functions) 
(natural exponential functions) 

(logarithmic functions) 
(natural logarithmic function). 


(summation rule) 
(product rule) 

(quotient rule) 

(chain rule) 

(logarithmic differentiation) 


, if / is one-to-one and onto. 


(differentiation of inverse functions). 


The methods of differential ealeulus just introdueed shall now be employed to describe the loeal 
change behaviour of a few simple examples of funetions in economic theory, and also to deter¬ 
mine their loeal extremal values. The following seetion provides an overview of sueh frequently 
oeeurring economic functions. 











1.3. COMMON FUNCTIONS IN ECONOMIC THEORY 

7.3 Common functions in economic theory 

1. total cost function K{x) >0 

argument: level of physieal output a: > 0 (dim: units) 

2. marginal cost function K'{x) > 0 

argument: level of physieal output x > 0 (dim: units) 

3. average cost function /f(x)/x > 0 

argument: level of physieal output x > 0 (dim: units) 

4. unit price function p{x) > 0 

argument: level of physieal output x > 0 (dim: units) 

5. total revenue function E{x) := xp{x) > 0 

argument: level of physieal output x > 0 (dim: units) 

6. marginal revenue function E'{x) = xp'{x) + p{x) 

argument: level of physieal output x > 0 (dim: units) 

7. profit function G{x) := E{x) — K{x) 
argument: level of physieal output x > 0 (dim: units) 
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(dim: CU) 

(dim: CU/unit) 
(dim: CU/unit) 
(dim: CU/unit) 
(dim: CU) 
(dim: CU/unit) 

(dim: CU) 


8. marginal profit function G"(x) := E'{x) — K\x) = xp'(x) + p(x) — K\x) (dim: 

CU/unit) 

argument: level of physieal output x > 0 (dim: units) 

9. utility function [/(x) (dim: case dependent) 

argument: material wealth, opportunity, action x (dim: case dependent) 

The fundamental concept of a utility function as a means to capture in quantitative terms 
the psychological value (happiness) assigned by an economic agent to a certain amount 
of money, or to owning a specific good, was introduced to economic theory in 1738 by 
the Swiss mathematician and physicist Daniel Bernoulli FRS (1700-1782); cf. Bernoulli 
(1738) JU. The utility function is part of the folklore of the theory, and often taken to be 
a piecewise differentiable, right-handedly curved (concave) function, i.e., U"{x) < 0, on 
the grounds of the assumption of diminishing marginal utility (happiness) with increasing 
material wealth. 


10. economic efficiency W{x) := E{x)/K{x) > 0 (dim: 1) 

argument: level of physical output x > 0 (dim: units) 

11. demand function N{p) > 0, monotonously decreasing (dim: units) 

argument: unit price p, (0 < p < Pmax) (dim: CU/unit) 

12. supply function A{p) > 0, monotonously increasing 
argument: unit price p, (pmin < P) (dim: CU/units). 


(dim: units) 
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A particularly prominent example of a real-valued eeonomie funetion of one real variable 
eonstitutes the psychological value function, devised by the Israeli-US-Ameriean experi¬ 
mental psyehologists Daniel Kahneman and Amos Tversky (1937-1996) in the eontext of 
their Prospect Theory (a pillar of Behavioural Economics), which was later awarded a 
Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel in 2002 (cf. Kahne¬ 
man and Tversky (1979) IfT^ p 279], and Kahneman (2011) IfTTl p 282f]). A possible representa¬ 
tion of this function is given by the piecewise description 

{ a log^o (1 + ^ ^ ^>0 

, (7.14) 

—2alogio (1 — a;) for x G M<o 

with parameter a G M>o. Overcoming a conceptual problem of Bernoulli’s utility function, here, 
in contrast, the argument x quantifies a change in wealth (or welfare) with respect to some given 
reference point (rather than a specific value of wealth itself). 


7.4 Curve sketching 

Before we turn to discuss applications of differential calculus to simple quantitative problems 
in economic theory, we briefly summarize the main steps of curve sketching for a real-valued 
function of one real variable, also referred to as analysis of the properties of differentiability of a 
real-valued function. 

1. domain: D{f) = {x G M|/(x) is regular} 

2. symmetries: for all x G D{f), is 

(i) /(-x) = /(x), i.e., is / even, or 

(ii) /(-x) = -fix), i.e., is / odd, or 

(iii) /(—x) 7 ^ /(x) 7 ^ —f{x), i.e., / exhibits no symmetries? 

3. roots: identify all xn G Dif) that satisfy the condition /(x) = 0. 

4. local extremal values: 

(i) local minima of / exist at all x^ G /?(/), for which the 

necessary condition /'(x) = 0, and the 

! 

sufficient condition /"(x) > 0 are satisfied simultaneously. 

(ii) local maxima of / exist at all x^ G D{f), for which the 

necessary condition /'(x) = 0, and the 

! 

sufficient condition /"(x) < 0 are satisfied simultaneously. 
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5. points of inflection: find all xw £ D{f), for which the 
necessary condition f"{x) = 0, and the 

I 

sufficient condition f"'{x) ^ 0 are satisfied simultaneously. 

6. monotonous behaviour: 

(i) / is monotonously increasing for all x G with j'{x) > 0 

(ii) / is monotonously decreasing for all x G with j'[x) < 0 

7. local curvature: 

(i) / behaves left-handedly curved for x G D{f) with f"{x) > 0 

(ii) / behaves right-handedly curved for x G D{f) with f"{x) < 0 

8. asymptotic behaviour: 

asymptotes to / are constituted by 

(i) straight lines y = ax + h with the property limj._j.+oo[/(x) — ax — h] = 0 or 
\im^^_^[f{x) - ax -b] = 0 

(ii) straight lines x = Xq at poles Xq ^ D{f) 

9. range: W{f) = {y e M|?/ = /(x)}. 


7.5 Analytic investigations of economic functions 


7.5.1 Total cost functions according to Turgot and von Thiinen 

According to the law of diminishing returns, which was introduced to economic theory by the 
French economist and statesman Anne Robert Jacques Turgot (1727-1781) and also by the Ger¬ 
man economist Johann Heinrich von Thiinen (1783-1850), it is meaningful to model non-negative 
total cost functions K (x) (in CU) relating to typical production processes, with argument level 
of physical output x > 0 units, as a mathematical mapping in terms of a special polynomial of 
degree 3 [cf. Eq. (17.51) 1. which is given by 


K (x) = a^x^ + a 2 X^ + Uix -f Oq 

=Kv{x) =Kf 


with as, ai > 0, 02 < 0, oq > 0, a| — SasOi < 0 . 


(7.15) 


The model thus contains a total of four free parameters. It is the outcome of a systematic regres¬ 
sion analysis of agricultural quantitative-empirical data with the aim to describe an inherently 
non-linear functional relationship between a few economic variables. As such, the functional 
relationship for K (x) expressed in Eq. (17.151) was derived from a practical consideration. It is a 
reflection of the following observed features: 
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(i) for levels of physieal output x > 0 units, the total eosts relating to typieal production 
processes exhibit strictly monotonously increasing behaviour; thus 

(ii) for the total costs there do not exist neither roots nor local extremal values however, 

(iii) the total costs display exactly one point of inflection. 

The continuous curve for K (x) resulting from these considerations exhibits the characteristic shape 
of an inverted capital letter “S”: beginning at a positive value corresponding to fixed costs, the total 
costs first increase degressively up to a point of inflection, whereafter they continue to increase, 
but in a progressive fashion. 

In broad terms, the functional expression given in Eq. (17.151) to model totals costs in dependence 
of the level of physical output is the sum of two contributions, the variable costs K^{x) and the 
fixed costs Kf = oq, viz. 

K{x) = K,{x)+Kj . (7.16) 

In economic theory, it is commonplace to partition total cost functions in the diminishing returns 
picture mio four phases, the boundaries of which are designated by special values of the level of 
physical output of a production process: 

• phase I (interval 0 units < x < xwY 

the total costs K{x) possess at a level of physical output xw = — 02 /( 303 ) > 0 units a 
point of inflection. For values of x smaller than xw, one obtains K"{x) < 0 CU/unit^, i.e., 
K (x) increases in a degressive fashion. For values of x larger than xw, the opposite applies, 
K"{x) > 0 CU/unit^, i.e., K{x) increases in a progressive fashion. The marginal costs, 
given by 


K'{x) = Soax^ + 2 a 2 X + Oi > 0 CU/unit for all x > 0 units , (7.17) 

attain a minimum at the same level of physical output, xw = — 02 /( 303 ). 

• phase II (interval xw < x <Xgf)\ 

the variable average costs 

KJx) 

— - -= oax^ + a 2 X + oi , x > 0 units (7.18) 

X 

become minimal at a level of physical output Xg^ = — 02 /( 203 ) > 0 units. At this value 
of X, equality of variable average costs and marginal costs applies, i.e., 

= K'{x) , (7.19) 

X 

^The last condition in Eq. (17.15b ensures a first derivative of K{x) that does not possess any roots; cf. the case of 
a quadratic algebraic equation 0 = ax^ -\-hx -\- c, with discriminant — 4oc < 0. 
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which follows by the quotient rule of differentiation from the neeessary eondition for an 
extremum of the variable average eosts, 

I / K^{x)V _ {K{x) - Kf)' ■ X - K^{x) ■ 1 
V a; J “ a;2 

and the fact that K'j: = 0 CU/unit. Taking care of the equality (17.191) . one finds for the tangent 
to K{x) in the point {xg^,K{xg^)) the equation [of. Eq. (17.131) 1 

T{x) = K{xg^)+K\xg^){x-Xg^) = K^{xg^)EKf + ^ {x-Xg^) = X . 

Xg-^ Xg^ 

Its interoept with the iT-axis is at iT/. 


• phase III (interval Xg^ < x < Xg^): 

The average costs 


Kix) 2 Oo 

-= Oj^fX T (X2X T cii T — , 

X X 


X > 0 units 


(7.20) 


attain a minimum at a level of physioal output Xg^ > 0 units, the defining equation of whioh 

is given by 0 CU = 2a3Xg^ + a 2 x‘^g^ — Oo- At this value of x, equality of average costs and 
marginal costs obtains, viz. 

^ = K'{x) , (7.21) 

X 

whioh follows by the quotient rule of differentiation from the neeessary eondition for an 
extremum of the average eosts. 


fK{x)\ K\x)-x-K{x)- 

\ X ) x2 


Since a quotient oan be zero only when its numerator vanishes (and its denominator remains 
non-zero), one finds from re-arranging the numerator expression equated to zero the property 


K'{x) 

K{x)/x "" K{x) 


K'{x) 

= X ' = 1 


for X = X, 


92 


(7.22) 


The eorresponding extremal value pair {xg^ , K {xg^ )) is referred to in economic theory as the 
minimum efficient scale (MES). From a business eeonomies perspeetive, at a level of phys¬ 
ioal output the (oompared to our remarks in the Introduotion inverted) ratio “INPUT 

K{x) 


over OUTPUT,” i.e.. 


X 


-, beeomes most favourable. By respecting the property (17.211) . 


the equation for the tangent to K (x) in this point [ef. Eq. (17.131) 1 beeomes 


T(x) = K{Xg^) + K'{Xg^){X - Xgfj = K {X gfj T 


K{x 


92 J 


X 


(X-XgJ = 


K{Xg^) 


X 


92 


X 


92 


Its interoept with the iT-axis is thus at 0 CU. 
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• phase IV (half-interval x > Xg^)\ 

In this phase K'{x)/K{x)/x > 1 obtains; the eosts assoeiated with the produetion of one 
additional unit of a good, approximately the marginal eosts K'{x), now exeeed the aver¬ 
age eosts, K{x)/x. This situation is eonsidered unfavourable from a business eeonomies 
perspeetive. 

7.5.2 Profit functions in the diminishing returns picture 

In this section, we confine our considerations, for reasons of simplicity, to a market sitution with 
only a single supplier of a good in demand. The price policy that this single supplier may thus 
inact defines a state of monopoly. Moreover, in addition we want to assume that for the market 
situation considered economic equilibrium obtains. This manifests itself in equality of supply 
and demand, viz. 

x{p) = N{p) , (7.23) 

wherein x denotes a non-negative supply function (in units) (which is synonymous with the sup¬ 
plier’s level of physical output) and N a non-negative demand function (in units), both of which 
are taken to depend on the positive unit price p (in CU/unit) of the good in question. The supply 
function, and with it the unit price, can, of course, be prescribed by the monopolistic supplier in 
an arbitrary fashion. In a specific quantitative economic model, for instance, the demand func¬ 
tion x{p) (recall that by Eq. (17.231) x{p) = N{p) obtains) could be assumed to be either a linear or 
a quadratic function of p. In any case, in order for x{p) to realistically describe an actual demand- 
unit price relationship, it should be chosen as a strictly monotonously decreasing function, and as 
such it is invertible. The non-negative demand function x{p) features two characteristic points, 
signified by its intercepts with the x- and the p-axes. The prohibitive price Pproh is to be deter¬ 
mined from the condition x(pproh) = 0 units; therefore, it constitutes a root of x{p). The saturation 
quantity Xsat, on the other hand, is defined by Xsat := a:(0 CU/unit). 

The inverse function associated with the strictly monotonously decreasing non-negative demand 
function x{p), the unit price function p{x) (in CU/unit), is likewise strictly monotonously 
decreasing. Via p{x), one calculates, in dependence on a known amount x of units sup¬ 
plied/demanded (i.e., sold), the total revenue (in CU) made by a monopolist according to (cf. 
Sec. 1731) 

E{x) = xp{x) . (7.24) 

Under the assumption that the non-negative total costs K (x) (in CU) underlying the production 
process of the good in demand can be modelled according to the diminishing returns picture of 
Turgot and von Thiinen, the profit function (in CU) of the monopolist in dependence on the level 
of physical output takes the form 


G{x) = E{x) - K{x) 


X p{x) — [a^x^ -\- a 2 x‘^ -t- aix 


+ 


total revenue 


total costs 


(7.25) 
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The first two derivatives of G{x) with respeet to its argument x are given by 

G'{x) = E'{x) — K'{x) = xp\x) + p{x) — + 2a2X + oi] (V.26) 

G'\x) = E'\x) — K'\x) = xp'\x) + 2p'{x) — [Qa^x + 202 ] . (7.27) 

Employing the principles of curve sketching set out in Sec. l7.4[ the following characteristic values 
of G{x) can thus be identified: 

• break-even point 

0:5 > 0 units, as the unique solution to the conditions 

G{x) = 0 CU (necessary condition) (7.28) 

and 

G'(a;) > 0 CU/unit (sufficient condition), (7.29) 

• end of the profitable zone 

xg > 0 units, as the unique solution to the conditions 

G{x) = 0 CU (necessary condition) (7.30) 

and 

G'(a;) < 0 CU/unit (sufficient condition), (7.31) 

• maximum profit 

xm > 0 CU, as the unique solution to the conditions 

G'(x) = 0 CU/unit (necessary condition) (7.32) 

and 

G"(a;) < 0 CU/unit^ (sufficient condition). (7.33) 

At this point, we like to draw the reader’s attention to a special geometric property of the quan¬ 
titative model for profit that we just have outlined: at maximum profit, the total revenue func¬ 

tion E{x) and the total cost function K{x) always possess parallel tangents. This is due to the 
fact that by the necessary condition for an extremum to exist, one finds that 

0 CU/unit = G'{x) = E\x) - K\x) ^ E\x) = K\x) . (7.34) 

GDC: Roots and local maxima resp. minima can be easily determined for a given stored function 
in mode CALC by employing the interactive routines zero and maximum resp. minimum. 

To conclude these considerations, we briefly turn to elucidate the technical term Cournot’s point, 
which frequently arises in quantitative discussions in economic theory; this is named after the 
French mathematician and economist Antoine-Augustin Cournot (1801-1877). Cournot’s point 
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simply labels the profit-optimal combination of the level of physical output and the associated unit 
price, {xm,p{xm)), for the unit price function p{x) of a good in a monopolistic market situation. 
Note that for this specific combination of optimal values the Amoroso-Robinson formula applies, 
which was developed by the Italian mathematician and economist Luigi Amoroso (1886-1965) 
and the British economist Joan Violet Robinson (1903-1983), This states that 


p{xm) 


K’{xm) 

1 + Ep{xm) 


(7.35) 


with K'{xm) the value of the marginal costs at xm, and ep{xM) the value of the elasticity of the 
unit price function at xm (see the following Sec. 17.61) . Starting from the defining equation of the 
total revenue E{x) = xp{x), the Amoroso-Robinson formula is derived by evaluating the first 
derivative of E{x) at xm, so 


E\xm) = p{xm) + xmp'{xm) = p{xm) 


1 + Xm 


p'{xm) 

p{xm) 


Sec. l7.6l 

= p{xm) [1 + £p{xm)] , 


and then re-arranging to solve for using the fact that E\xm) = K\xm)- 

Remark: In a market situation where perfect competition applies, one assumes that the unit 
price function has settled to a constant value p{x) = p = constant > 0 CU/unit (and, hence, 
p'{x) = 0 CU/unit^ obtains). 


7.5.3 Extremal values of rational economic functions 


Now we want to address the determination of extremal values of economic functions that constitute 
ratios in the sense of the construction 

OUTPUT 
INPUT ’ 

a topic raised in the Introduction. 

Let us consider two examples for determining local maxima of ratios of this kind. 


(i) We begin with the average profit in dependence on the level of physical output x > 0 units. 


G{x) 


X 


(7.36) 


The conditions that determine a local maximum are [G{x)/x\' = 0 CU/unit^ and 

[G{x)/x\" < 0 CU/unit^. Respecting the quotient rule of differentiation (cf. Sec. 17.21) . 
the first condition yields 


G'{x)x — G{x) 


x^ 


= 0 GE/ME^ . 


(7.37) 


Since a quotient can only be zero when its numerator vanishes while its denominator remains 
non-zero, it immediately follows that 


G\x)x - G{x) = 0 CU 


G\x) 

1 

G[x) 


(7.38) 
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The task at hand now is to find a (unique) value of the level of physieal output x whieh 
satisfies this last eondition, and for whieh the seeond derivative of the average profit beeomes 
negative. 


(ii) To eompare the performanee of two eompanies over a given period of time in a meaningful 
way, it is recommended to adhere only to measures that are dimensionless ratios, and so 
independent of scale. An example of such a dimensionless ratio is the measure referred to 

as economic efficiency, 

Eix) 

^ , (V.39) 

K[x) 

which expresses the total revenue (in CU) of a company for a given period as a multiple of 
the total costs (in CU) it had to endure during this period, both as functions of the level of 
physical output. In analogy to our discussion in (i), the conditions for the existence of a 

local maximum amount to [E{x)/K{x)]' = 0 x 1/unit and [E{x)/K{x)]" < 0 x 1/unit^. 
By the quotient rule of differentiation (see Sec. l7.2L the first condition leads to 


E\x)K{x) - E{x)K\x) 
A2(x) 


0 X 1/unit, 


(7.40) 


i.e., ior K{x) > 0 CU, 


E'{x)K{x) - E{x)K'{x) = 0 CUVunit. 


(7.41) 


By re-arranging and multiplication with x > 0 unit, this can be cast into the particular form 


E\x) K\x) 
E{x) ~ ^ K{x) ■ 


(7.42) 


The reason for this special kind of representation of the necessary condition for a local 
maximum to exist [and also for Eq. (17.381) 1 will be clarified in the subsequent section. Again, 
a value of the level of physical output which satisfies Eq. (17.421) must in addition lead to 
a negative second derivative of the economic efficiency in order to satisfy the sufficient 
condition for a local maximum to exist. 


7.6 Elasticities 

Einally, we pick up once more the discussion on quantifying the local variability of differentiable 
real-valued functions of one real variable, f : D CK ^ W CM, though from a slightly different 
perspective. Eor reasons to be elucidated shortly, we confine ourselves to considerations of regimes 
of / with positive values of the argument x and also positive values y = f{x) > 0 of the function 
itself. 

As before in Sec. 17.2[ we want to assume a small change of the value of the argument x and 
evaluate its resultant effect on the value y = f{x). This yields 

X + Ax y = f{x) y + Ay = f{x + Ax) . (7.43) 


X 
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We remark in passing that relative changes of non-negative quantities are defined by the quotient 

new value — old value 
old value 

under the prerequisite that “old value > 0” applies. It follows from this specific construction that 
the minimum value a relative change can possibly attain amounts to “—1” (corresponding to a 
decrease of the “old value” by 100%). 

Related to this consideration we identify the following terms: 


• a prescribed absolute change of the independent variable x: Ax , 

• the resultant absolute change of the function /: Ay = f{x + Ax) — f{x) , 

Ax 

• the associated relative change of the independent variable x: — , 

X 


• the associated resultant relative change of the function /: 


^ + Ax) - f{x) 

y f{x) 


/\ X 

Now let us compare the order-of-magnitudes of the two relative changes just envisaged, — and 

Ay 
y 


X 


. This is realised by considering the value of their quotient, “resultant relative change of / 


divided by the prescribed relative change of x’”: 

Ay f{x + Ax) - f{x) 


y 

Ax 

X 


fix) 


Ax 

X 


Since we assumed / to be differentiable, it is possible to investigate the behaviour of this quotient 

Ax 

of relative changes in the limit of increasingly smaller prescribed relative changes - —)■ 0 

Ax —)■ 0 near some a; > 0. One thus defines: 


X 


Def.: For a differentiable real-valued function / of one real variable x, the dimensionless (i.e., 
units-independent) quantity 


ef{x) := lim 


Ay 

y 


Ax 


lim ■ 

Aa;—>-0 


f{x + Ax) - f{x) 

fix) 

Ax 


= X 


fix) 

fix) 


(7.44) 


X 


X 


is referred to as the elasticity of the function / at position x. 

The elasticity of / quantifies its resultant relative change in response to a prescribed infinitesimally 
small relative change of its argument x, starting from some positive initial value x > 0. As such it 
constitutes a measure for the relative local rate of change of a function / in a point (x,/(x)). In 
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economic theory, in particular, one adheres to the following interpretation of the elastieity ef{x): 
when the postive argument x of some positive differentiable real-valued funetion / is inereased by 
1 %, then in eonsequenee / will ehange approximately by ef{x) x 1 %. 

In the scientific literature one often finds the elasticity of a positive differentiable funetion / of a 
positive argument x expressed in terms of logarithmie differentiation. That is, 

Sfix) := for X > 0 and f{x) > 0 , 

^ dln(x) ^ ^ 

since by the chain rule of differentiation it holds that 


d/(a^) d/(x) 

^ fix) _ ~dfr _ fi^) 

dln(x) dx /(x) /(x) 

X 

The logarithmie representation of the elastieity of a differentiable funetion / immediately explains 
why, at the beginning, we confined our considerations to positive differentiable functions of posi¬ 
tive arguments onlyO A brief look at the list of standard eeonomie funetions provided in Sec. 17.31 
reveals that most of these (though not all) are positive funetions of non-negative arguments. 

For the elementary elasses of real-valued funetions of one real variable diseussed in See. 17. H one 
finds: 

Standard elasticities 

1. /(x) = x*^ for n e N and x G M>o ef{x) = n (natural power-law functions) 

2. /(x) = x" for a G M and x G M>o ^ ef{x) = a (general power-law functions) 

3. /(x) = a* for a G M>o\{l} and x G M>o ^f{x) = ln(a)x (exponential functions) 

4. /(x) = for a G M and x G M>o ^ ef(x) = ax (natural exponential functions) 

5. /(x) = logQ(x) for a G M>o\{l} and x G M>o 

^ Sfix) = , , (logarithmic functions) 

ln(a)log„(x) 


6 . /(x) = ln(x) forx G M>o 



(natural logarithmic function). 


In view of these results, we would like to emphasise the faet that for the entire family of gen¬ 
eral power-law functions the elastieity £/(x) has a constant value, independent of the value of 
the argument x. It is this very property whieh elassifies general power-law functions as scale- 
invariant. When scale-invariance obtains, dimensionless ratios, i.e., quotients of variables of the 

^To extend the regime of applicability of the measure e /, one may consider working in terms of absolute values |a;| 
and |/(a;) |. Then one has to distinguish between four cases, which need to be looked at separately: (i) a; > 0, f{x) > 0, 
(ii) a: < 0, /(x) > 0, (iii) x < 0, /(x) < 0 and (iv) x > 0, /(x) < 0. 
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same physical dimension, reduce to constants. In this context, we would like to remark that scale- 
invariant (fractal) power-law functions of the form f{x) = Kx^^, with K > t) and ot G M<o\{—1}, 
are frequently employed in Economics and the Social Sciences for modelling uncertainty of eco¬ 
nomic agents in decision-making processes, or for describing probability distributions of rare 
event phenomena; see, e.g., Taleb (2007) p 326ff] or Gleick (1987) [[TSl Chs. 5 and 6]. This 
is due, in part, to the curious property that for certain values of the exponent a general power-law 
probability distributions attain unbounded variance; cf. Ref. IfT^ Sec. 8.9]. 


Practical applications in economic theory of the concept of an elasticty as a measure of relative 
change of a differentiable real-valued function / of one real variable x are generally based on the 
following linear (!) approximation: beginning at Xq > 0, for small prescribed percentage changes 

/\t. 

of the argument x in the interval 0 % <-< 5 %, the resultant percentage changes of / amount 

Xo 


approximately to 


(percentage change of /) (elasticity of / at xq) x (percentage change of x) , (7.45) 


or, in terms of a mathematical formula, to 


/(xq -f Nx) - f(xo) 
f{xo) 


Efixo) 


Ax 

Xo 


(7.46) 


We now draw the reader’s attention to a special kind of terminology developed in economic theory 
to describe the relative local change behaviour of economic functions in qualitative terms. For 
X G D(f), the relative local change behaviour of a function / is called 

• inelastic, whenever \sf{x)\ < 1, 

• unit elastic, when |£/(x)| = 1, and 

• elastic, whenever I £/(x) I > 1. 

For example, a total cost function K (x) in the diminishing returns picture exhibits unit elastic 
behaviour at the minimum efficient scale where, by Eq. (I7.22L exixgO = 1. Also, at 

the local maximum of an average profit function G{x)/x, the property Scix) = 1 applies; cf. 

Eq. (17381) . 

Next, we review the computational rules one needs to adhere to when calculating elasticities for 
combinations of two real-valued functions of one real variable in the sense of Sec. 17.1.61 

Computational rules for elasticities 

If / and g are differentiable real-valued functions of one real variable, with elasticities Sf and Sg, 
it holds that: 

1. product / ■ g: = £f{x) + eg{x), 

f 

2. quotient -,gy^0: £f/g{x) = £/(x) - £g(x), 

3. concatenation fog: £fog{x) = ef{g{x)) ■ £g{x). 
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4. inverse function / 




1 

^fiy) 


y=f Ha;) 


To end this chapter, we remark that for a positive differentiable real-valued function / of one 
positive real variable x, a second elasticity may be defined according to 




d X 

^ dx f{x) 


d/(a^) 

dx 


(7.47) 


Of course, by analogy this procedure may be generalised to higher derivatives of / still. 
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Chapter 8 


Integral calculus of real-valued functions of 
one real variable 


In the final chapter of these lecture notes we give a brief overview of the main definitions and 
laws of the integral calculus of real-valued functions of one variable. Subsequently we consider 
a simple application of this tool in economic theory. 

8.1 Indefinite integrals 

Def.: Let / be a continuous real-valued function of one real variable and F a differentiable real¬ 
valued function of the same real variable, with D{f) = D{F). Given that / and F are related 
according to 


F'{x) = f{x) for all x e D{f) , 


( 8 . 1 ) 


then F is referred to as a primitive off. 

Remark: The primitive of a given continuous real-valued function / cannot be unique. By the 
rules of differentiation discussed in Sec. 17.21 besides F also F + c, with c G M a real-valued 
constant, constitutes a primitive of / since (c)' = 0. 

Def.: If F is a primitive of a continuous real-valued function / of one real variable, then 



c = constant e M, •with F'(x) = f (x) (8.2) 


defines the indefinite integral of the function /. The following names are used to refer to the 
different ingredients in this expression: 


X — the integration variable, 

f{x) — the integrand, 

da; — the differential, and, lastly. 
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• c — the constant of integration. 

For the elementary, continuous real-valued functions of one variable introduced in Sec. 17.11 the 
following rules of indefinite integration apply: 

Rules of indefinite integration 

1. J a dx = ax + c with a = constant G M 

„2 


2 . 


X 


xdx = "^ + c 


(constants) 
(linear functions) 


3. x^dx = 


X 


,n+l 


n-\-l 


4. 


dx = 


X 


.a-l-l 


5. / dx = 


Q: -j- 1 
a* 

ln(a) 


-f c forn G N 

-f c for a G M\{—1} and x G 
-f c for a G M>o\{l} 


6. / 6°“^ dx = -h c for a G M\{0} 


(natural power-law functions) 
&>o (general power-law functions) 
(exponential functions) 
(natural exponential functions) 


7. J X ^ dx = In |x| -f c for x G M\{0}. 

Special methods of integration need to be employed when the integrand consists of a concatanation 
of elementary real-valued functions. Here we provide a list with the main tools for this purpose. 
For differentiable real-valued functions / and g, it holds that 

1. /(a/(x) ± f3g{x)) dx = a f /(x) dx ± /3 f g(x) dx 

with a, I3 = constant G M (summation rule) 

2. J f{x)g'{x) dx = f{x)g{x) — J f'{x)g{x) dx (integration by parts) 

u=g(x) and du=g'{x)dx 

3. J f{g{x))g'{x)dx J f{u) du = F{g{x))c (substitution method) 


f'{x) 

fix) 


dx = In |/(x)| -f c for /(x) ^ 0 


(logarithmic integration). 


8.2 Definite integrals 


Def.: Let / be a real-valued function of one variable which is continuous on an interval [a, b] C 
D{f), and let F be a primitive of /. Then the expression 


(8.3) 



defines the definite integral of / in the limits of integration a and b. 

For definite integrals the following general rules apply: 
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1 . 

2 . 

3. 


f{x) dx 

I 

f{x) dx 

I 

/(x) dx 


0 


(identical limits of integration) 



/(x) dx 



/(x) dx + 



(interchange of limits of integration) 
x) dx for c e [a, b] (split of integration interval). 


Remark; The main qualitative differenee between an (i) indefinite integral and a (ii) definite in¬ 
tegral of a eontinuous real-valued function of one variable reveals intself in the different kinds of 
outcome: while (i) yields as a result a real-valued (primitive)/Mnchon, (ii) simply yields a single 
real number. 


GDC: For a stored real-valued function, the evaluation of a definite integral can be performed in 
mode CALC with the pre-programmed function Jf (x) dx. The corresponding limits of integration 
need to be specified interactively. 


As indicated in Sec. 17.61 the scale-invariant power-law functions /(x) = 
X e M>o play a special role in practical applications. For x G [a, h] C M>o 
that 



X 


a+1 


Of -|- 1 


x=b 


Of 1 




Q+l 


a“+i) 


x“ for a G M and 
and a ^ —1 it holds 

(8.4) 


Problematic in this context can be considerations of taking limits of the form a —?• 0 resp. b ^ oo, 
since for either of the two cases 


(i) case a < —1: 


(ii) case a > — 1: 


lim / x" dx ^ oo , 
a^>0 I 


lim / x" dx —)■ cxo , 

6—>-oo I^ 


one ends up with divergent mathematical expressions. 


(8.5) 


( 8 . 6 ) 


8.3 Applications in economic theory 

The starting point shall be a simple market situation for a single product. For this product, on the 
one-hand side, there be a demand function N (p) (in units) which is monotonously decreasing on 
the price interval [p„, po]; the limit values p„ and po denote the minimum and maximum prices per 
unit (in CU/u) acceptable for the product. On the other hand, the market situation be described by 
a supply function A{p) (in units) which is monotonously increasing on [p^J,Po]■ 

The equilibrium unit price pm (in CU/unit) for this product is defined by assuming a state of 
economic equilibrium of the market, quantitatively expressed by the condition 


A{pm) = N{pm) ■ 


( 8 . 7 ) 
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Geometrically, this condition defines common points of intersection for the functions A{p) and 
N{p) (when they exist). 

GDC: Common points of intersection for stored functions / and g can be easily determined inter¬ 
actively in mode CALC employing the routine intersect. 

In a drastically simplified fashion, we now turn to compute the revenue made on the market by the 
suppliers of a new product for either of three possible strategies of market entry. 

1. Strategy 1: The revenue obtained by the suppliers when the new product is being sold 
straight at the equilibrium unit price pm, in an amount N{pm), is simply given by 

Ui = U{pm) = PmN{pm) (inCU) . (8.8) 


2. Strategy 2: Some consumers would be willing to purchase the product intially also at a unit 
price which is higher than pm- If, hence, the suppliers decide to offer the product initially 
at a unit price Po > Pm, and then, in order to generate further demand, to continuousljQ (!) 
reduce the unit price to the lower pm, the revenue obtained yields the larger value 

fPo 

U2 = U{pm)+ N{p)dp. (8.9) 

JPM 

Since the amount of money 

fPo 

K:= N{p)dp (inCU) (8.10) 

•I Pm 

is (theoretically) safed by the consumers when the product is introduced to the market ac¬ 
cording to strategy 1, this amount is referred to in the economic literature as consumer 
surplus. 


3. Strategy 3: Some suppliers would be willing to introduce the product to the market initially 
at a unit price which is lower than pm- If, hence, the su|mliers decide to offer the product 
initially at a unit price pu < Pm, and then to continuouslja (!) raise it to the higher pm, the 
revenue obtaines amounts to the smaller value 


= U{pm) 



A{j)) dp . 


Since the suppliers (theoretically) earn the extra amount 


P := 



A{p) dp 


(in CU) 


( 8 . 11 ) 


( 8 . 12 ) 


when the product is introduced to the market according to strategy 1, this amount is referred 
to in the economic literature as producer surplus. 


'This is a strong mathematical assumption aimed at facilitating the actual calculation to follow. 
^See previous footnote. 




Appendix A 

Glossary of technical terms (GB - D) 


A 

absolute change: absolute Anderung 
absolute value: Betrag 
account balance: Kontostand 
addition: Addition 

analysis: Analysis, Untersuchung auf Differenzierbarkeitseigenschaften 

arithmetical mean: arithmetischer Mittelwert 

arithmetical sequence: arithmetische Zahlenfolge 

arithmetical series: arithmetische Reihe 

augmented coefficient matrix: erweiterte Koeffizientenmatrix 

average costs: Stiickkosten 

average profit: Durchschnittsgewinn, Gewinn pro Stiick 

B 

backward substitution: riickwertige Substitution 
balance equation: Bilanzgleichung 
basis: Basis 

basis solution: Basisldsung 
basis variable: Basisvariable 
Behavioural Economics: Verhaltensokonomik 
boundary condition: Randbedingung 
break-even point: Gewinnschwelle 

C 

chain rule: Kettenregel 

characteristic equation: charakteristische Gleichung 

coefficient matrix: Koeffizientenmatrix 

column: Spalte 

column vector: Spaltenvektor 

component: Komponente 

compound interest: Zinseszins 

concatenation: Verschachtelung, Verkniipfung 

conservation law: Erhaltungssatz 
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constant of integration: Integrationskonstante 
constraint: Zwangsbedingung 
cost function: Kostenfunktion 
consumer surplus: Konsumentenrente 
continuity: Stetigkeit 
contract period: Laufzeit 
Cournot’s point: Cournotscher Punkt 
curve sketching: Kurvendiskussion 

D 

decision-making: Entscheidungsfindung 

declining-balance depreciation method: geometrisch-degressive Abschreibung 

definite integral: bestimmtes Integral 

demand function: Nachfragefunktion 

dependent variable: abhangige Variable 

depreciation: Abschreibung 

depreciation factor: Abschreibungsfaktor 

derivative: Ableitung 

determinant: Determinante 

difference: Differenz 

difference quotient: Differenzenquotient 

differentiable: differenzierbar 

differential: Integrationsdifferenzial 

differential calculus: Differenzialrechnung 

dimension: Dimension 

direction of optimisation: Optimierungsrichtung 
divergent: divergent, unbeschr”ankt 
domain: Definitionsbereich 

E 

economic agent: Wirtschaftstreibende(r) (meistens ein homo oeconomicus) 

economic efficiency: Wirtschaftlichkeit 

economic equilibrium: okonomisches Gleichgewicht 

economic principle: okonomisches Prinzip 

economic theory: Wirtschaftstheorie 

Econophysics: Okonophysik 

eigenvalue: Eigenwert 

eigenvector: Eigenvektor 

elastic: elastisch 

elasticity: Elastizitat 

element: Element 

end of profitable zone: Gewinngrenze 
endogenous: endogen 
equilibrium price: Marktpreis 

equivalence transformation: Aquivalenztransformation 
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exogenous: exogen 

exponential function: Exponentialfunktion 

extrapolation: Extrapolation, iiber bekannten Gutigkeitsbereich hinaus verallgemeinern 

F 

feasible region: zulassiger Bereich 
final capital: Endkapital 
fixed costs: Eixkosten 
forecasting: Vorhersagen erstellen 
function: Eunktion 

G 

GauBian elimination: GauB’scher Algorithmus 
GDC: GTR, grafikfahiger Taschenrechner 
geometrical mean: geometrischer Mittelwert 
geometrical sequence: geometrische Zahlenfolge 
geometrical series: geometrische Reihe 
growth function: Wachstumsfunktion 

H 

I 

identity: Identitat 

image vector: Absolutgliedvektor 

indefinite integral: unbestimmtes Integral 

independent variable: unabhangige Variable 

inelastic: unelastisch 

initial capital: Anfangskapital 

installment: Ratenzahlung 

installment savings: Ratensparen 

integral calculus: Integralrechnung 

integrand: Integrand 

integration variable: Integrationsvariable 

interest factor: Aufzinsfaktor 

interest rate: ZinsfuB 

inverse function: Inversfunktion, Umkehrfunktion 
inverse matrix: inverse Matrix, Umkehrmatrix 
isoquant: Isoquante 

J 

K 

L 

law of diminishing returns: Ertragsgesetz 
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length: Lange 

level of physical output: Ausbringungsmenge 
limits of integration: Integrationsgrenzen 
linear combination: Linearkombination 
linearisation: Linearisierung 
linear programming: lineare Optimierung 
local rate of change: lokale Anderungsrate 
logarithmic function: Logarithmusfunktion 

M 

mapping: Abbildung 

marginal costs: Grenzkosten 

maximisation: Maximierung 

minimisation: Minimierung 

minimum efficient scale: Betriebsoptimum 

monetary value: Geldwert 

monopoly: Monopol 

mono tonicity: Monotonie 

mortgage loan: Darlehen 

N 

non-basis variable: Nichtbasisvariable 

non-negativity constraints: Nichtnegativitatsbedingungen 

non-linear functional relationship: nichtlineare Funktionalbeziehung 

O 

objective function: Zielfunktion 
one-to-one and onto: eineindeutig 
optimal solution: optimalen Losung 
optimal value: optimaler Wert 
optimisation: Optimierung 
order-of-magnitude: GroBenordnung 
orthogonal: orthogonal 
over-determined: iiberbestimmt 

P 

parallel displacement: Parallelverschiebung 
pension calculations: Rentenrechnung 
percentage rate: Prozentsatz 
perfect competition: to tale Konkurrenz 
period: Periode 

pivot column index: Pivotspaltenindex 
pivot element: Pivotelement 
pivot operation: Pivotschritt 
pivot row index: Pivotzeilenindex 
pole: Polstelle, Singularitat 
polynomial division: Polynomdivision 
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polynomial of degree n: Polynom vom Grad n 
power-law function: Potenzfunktion 
present value: Barwert 
primitive: Stammfunktion 

principal component analysis: Hauptkomponentenanalyse 

producer surplus: Produzentenrente 

product rule: Produktregel 

profit function: Gewinnfunktion 

prohibitive price: Pohibitivpreis 

Prospect Theory: Neue Erwartungstheorie 

psychological value function: psychologische Wertfunktion 

Q 

quadratic matrix: quadratische Matrix 

quotient: Quotient 

quotient rule: Quotientenregel 

R 

range: Wertespektrum 
rank: Rang 

rare event: seltenes Ereignis 
rational function: gebrochen rationale Eunktion 
real-valued function: reellwertige Eunktion 
reference period: Referenzzeitraum 
regression analysis: Regressionsanalyse 
regular: regul^ 

relative change: relative Anderung 
remaining debt: Restschuld 
remaining resources: Restkapazitaten 
remaining value: Restwert 
rescaling: Skalierung 
resources: Rohstoffe 

resource consumption matrix: Rohstoffverbrauchsmatrix 
restrictions: Restriktionen 
root: Nullstelle 
row: Reihe 

row vector: Zeilenvektor 

S 

saturation quantity: Sattigungsmenge 
scale: Skala, GroBenordnung 
scale-invariant: skaleninvariant 
simplex: Simplex, konvexer Polyeder 
simplex tableau: Simplextabelle 
singular: singula 
sink: Senke 
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slack variable: Schlupfvariable 
source: Quelle 

stationary: station^, konstant in der Zeit 

straight line depreciation method: lineare Abschreibung 

strictly monotonously decreasing: streng monoton fallend 

strictly monotonously increasing: streng monoton steigend 

summation rule: Summationsregel 

supply function: Angebotsfunktion 

T 

tangent: Tangente 

target space: Wertebereich 

technology matrix: Technologiematrix 

total demand matrix: Gesamtbedarfsmatrix 

total revenue: Ertrag 

transpose: Transponierte 

U 

uncertainty: Unsicherheit 
under-determined: unterbestimmt 
uniqueness: Eindeutigkeit 
unit elastic: proportional elastisch 
unit matrix: Einheitsmatrix 
unit price: Stiickpreis 
unit vector: Einheitsvektor 
utility function: Nutzenfunktion 

V 

value chain: Wertschopfungskette 
variability: Anderungsverhalten, Variabilitat 
variable average costs: variable Stiickkosten 
variable costs: variable Kosten 
variable vector: Variablenvektor 
vector: Vektor 

vector algebra: Vektoralgebra 

W 

well-determined: wohlbestimmt 

Z 

zero matrix: Nullmatrix 
zero vector: Nullvektor 



Bibliography 


[1] T Arens, F Hettlich, C Karpfinger, U Kockelkorn, K Lichtenegger and H Stachel (2008) 
Mathematik 1®* Edition (Heidelberg: Spektrum) ISBN-13: 9783827417589 

[2] A Asano (2013) An Introduction to Mathematics for Economics (Cambridge: Cambridge 
University Press) ISBN-13: 9780521189460 

[3] C Bauer, M Clausen, A Kerber and H Meier-Reinhold (2008) MathematikfUr Wirtschaftswis- 
senschaftler 5^'^ Edition (Stuttgart: Sehaffer-Poesehel) ISBN-13: 9783791027487 

[4] Bernoulli D (1738) Speeimen theoriae novae de mensura sortis Papers of the Imperial 
Academy of Sciences in St. Petersburg 

English translation: (1954) Exposition of a new theory on the measurement of risk 
Econometrica 22 23-36 

[5] J-P Bouehaud and M Potters (2003) Theory of Einancial Risk and Derivative Pricing 2“'^ 
Edition (Cambridge: Cambridge University Press) ISBN-13: 9780521741866 

[6] K Boseh (2003) Mathematik fur Wirtschaftswissenschaftler Edition (Miinehen: Olden- 

bourg) ISBN-10: 3486273094 

[7] IN Bronstein, K A Semedjajew, G Musiol and H Miihlig (2005) Taschenbuch der Mathematik 
6**^ Edition (Erankfurt (Main): Harri Deutseh) ISBN-10: 3817120060 

[8] G B Dantzig (1949) Programming of interdependent aetivities: II mathematieal model 
Econometrica 17 200-211 

[9] G B Dantzig (1955) Einear programming under uneertainty Management Science 1 197-206 

[10] E T Dowling (1990) Schaum’s Outline of Calculus for Business, Economics, and the Social 
Sciences (New York: McGraw-Hill) ISBN-13: 9780070176737 

[11] E T Dowling (2009) Schaum ’s Outline of Mathematical Methods for Business and Economics 
(New York: McGraw-Hill) ISBN-13: 9780071635325 

[12] H van Elst (2015) Eoundations of descriptive and inferential statistics Preprint 
arXiv:1302.2525v3 [stat.AP] 

[13] J Gleick (1987) Chaos — Making a New Science Edition 1998 (Eondon: Vintage) ISBN- 
13: 9780749386061 


81 



82 


BIBLIOGRAPHY 


[14] J F Hair jr, W C Black, B J Babin and R E Anderson (2010) Multivariate Data Analysis 7*^ 
Edition (Upper Saddle River (NJ): Pearson) ISBN-13: 9780135153093 

[15] E Herrmann (2003) Physik III: Thermodynamik — Skripten zur Experimentalphysik Univer- 
sitat Karlsruhe (TH) URL (eited on August 30, 2015): 
www.physikdidaktik.uni-karlsruhe.de/skripten/thermod.pdf 

[16] J Hiilsmann, W Gamerith, U Leopold-Wildburger and W Steindl (2005) Einfuhrung in die 
Wirtschaftsmathematik Edition (Berlin: Springer) ISBN-10: 3540244093 

[17] D Kahneman (2011) Thinking, East and Slow (London: Penguin) ISBN-13: 9780141033570 

[18] D Kahneman and A Tversky (1979) Prospeet theory: an analysis of deeision under risk 
Econometrica 47 263-292 

[19] S J Leon (2009) Linear Algebra with Applications 8*^ Edition (London: Pearson) ISBN-13: 
9780136009290 

[20] W W Leontief (1936) Quantitative input and output relations in the eeonomie systems of the 
United States The Review of Economics and Statistics 18 105-125 

[21] R Penrose (2004) The Road to Reality — A Complete Guide to the Laws of the Universe 1®* 
Edition (London: Jonathan Cape) ISBN-10: 0224044478 

[22] H Rinne (2008) Taschenbuch der Statistik 4*^ Edition (Erankfurt Main): Harri Deutseh) 
ISBN-13: 9783817118274 

[23] S Singh (1997) Eermat’s Last Theorem (London: Louth Estate) ISBN-13: 9781857025217 

[24] L Smith (2007) Chaos — A Very Short Introduction (Oxford: Oxford University Press) 
ISBN-13: 9780192853783 

[25] N N Taleb (2007) The Black Swan — The Impact of the Highly Improbable (London: Pen¬ 
guin) ISBN-13: 9780141034591 

[26] A Zeh-Marsehke (2010) Einanzmathematik — Zins und Zinseszins, Tilgungs- und Renten- 
rechnung, Investitionsrechnung, Abschreibungen private leeture notes 



